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ABSTRACT 


A  model  is  developed  which  provides  for  the  generation  of  the 
approximate  spectra  of  pulsed  signals  with  arbitrary  amplitude 
and  frequency  modulation.  Error  bounds  associated  with  the  model  are 
developed  providing  an  insight  into  performance  accuracy.  The  concept 
of  using  two  approximate  spectra  to  generate  a  more  exact  spectrum 
by  means  of  computationally  trivial  calculations  is  addressed. 

The  model  is  used  to  validate  the  feasibility  of  using  frequency  domain 
discriminants  in  the  signal  classification  process. 
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I.  INTRODUCTION 


The  classical  radar  pulse  consists  of  a  rectangular  pulse  modulating 
a  constant  carrier  frequency.  This  gives  rise  to  the  characteristic 
[sin(x)]/x  spectrum  centered  about  the  carrier  frequency. 

In  practice,  this  idealized  pulse  is  not  realized  since  pulse 
envelopes  tend  to  be  more  trapezoidal  than  rectangular  and  since  carrier 
frequencies  are  not,  de  facto,  constant.  Frequency  modulation  (FM), 
intentional  or  unintentional,  may  be  a  property  of  all  radars.  As  an 
example  of  this,  oscillators  may  experience  finite  turn-on  and  shut-down 
periods  during  which  a  transient  frequency  modulation  can  be  introduced 
onto  the  carrier.  Figure  1  illustrates  what  a  typical  instantaneous 
frequency  of  a  "constant"  carrier  radar  pulse  might  look  like. 

The  question  of  how  modulation  affects  the  frequency  spectrum  is 
one  of  interest  in  pulse  analysis.  Thus  it  is  desirable  to  model  the 
amplitude  spectrum  of  a  pulse  with  arbitrary  amplitude  and  frequency 
modulation  to  provide  a  vehicle  to  study  the  effects  of  modulation  on 
the  amplitude  spectrum. 

One  approach  to  this  general  modeling  problem  is  to  generate  a  digital 
representation  of  the  pulse  waveform  and  apply  a  Discrete  Fourier  Trans¬ 
form  (DFT)  algorithm  to  the  data  to  generate  the  frequency  domain 
representation.  This  approach  for  the  general  case  can  be  cumbersome 
in  generating  the  digital  signal,  satisfying  the  Nyquist  criteria,  and 
trading-off  between  frequency  resolution  and  transform  size. 

An  approach  to  addressing  the  amplitude  spectrum  of  frequency  modu¬ 
lation  waveforms  is  to  view  the  problem  in  terms  of  the  quasi-stationary 
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Figure  1  - 


tantaneous  Frequency  Function 
ociated  with  a  Pulsed  Signal 


r.Tfcii 


principle.  The  fraction  of  time  that  the  instantaneous  frequency  for  the 
waveform  spends  in  a  particular  frequency  interval  is  the  fraction  of 
power  that  belongs  to  the  corresponding  Fourier  frequency  interval. 

This  principle,  however,  is  only  validly  applicable  asymptotically  for 
large  modulation  indices.  [1]. 

Holway  and  Mullen  [1]  have  addressed  a  special  case  of  FM  spectra 
generation.  Using  a  Fourier  Series  approach  and  with  the  aid  of  Fresnel 
Integrals,  a  closed  form  expression  has  been  developed  for  the  line 
spectrum  for  periodic  frequency  modulated  signals  in  those  cases  when 
the  modulating  function  is  trapezoidal.  This  model  is  flexible  in  that 
the  form  of  the  trapezoidal  modulation  can  be  varied  from  square-wave 
modulation  at  one  extreme  to  saw-tooth  modulation  at  the  other. 

As  with  the  work  of  Holway  and  Mullen,  the  model  developed  in  this 
paper  uses  the  Fresnel  integral  as  one  of  its  keystones,  but  addresses 
a  broader  scope  of  modulated  signals.  The  Elemental  Cell  Model  which 
is  to  be  addressed  is  designed  to  generate  the  approximate  spectrum  of 
pulsed  signals  with  arbitrary  amplitude  and  frequency  modulation.  The 
approach  is  to  partition  the  amplitude  envelope  and  instantaneous 
frequency  function  associated  with  a  signal  into  an  arbitrary  number  of 
"elemental  cells"  of  arbitrary  length.  Across  each  elemental  cell  the 
amplitude  and  frequency  is  linearly  approximated  so  that  the  original 
envelope  and  instantaneous  frequency  function  have  both  been  formed  into 
piecewise  linear  continuous  functions.  Each  elemental  cell  has  associated 
with  it  a  characteristic  frequency  spectrum  which  can  be  represented  in 
closed  form.  With  phase  continuity  maintained  between  contiguous  cells, 
the  individual  elemental  cell  spectra  can  be  summed  to  provide  the 
approximated  spectrum  of  the  original  signal.  This  model  provides  a 


flexible,  easy  to  use  tool  for  spectrum  analysis. 


14 


II.  ELEMENTAL  CELL  MODEL 


A.  DEVELOPMENT  OP  THE  ELEMENTAL  CELL  MODEL  FOR  THE  UNIT  PULSE 


In  this  section,  the  elemental  cell  model  is  developed  £or  the  unit 
amplitude  pulse  waveform.  The  model  provides  a  flexible,  easily  applied 
method  of  generating  approximated  spectra  for  pulsed  signals  with  arbitrary 
frequency  modulation.  The  approach  of  the  model  is  to  decompose  pulse 
waveforms  into  sequences  of  elemental  (waveforms)  cells,  each  of  which 
has  a  characteristic  spectrum  in  closed  analytic  form. 

Consider  a  unit  pulse  with  an  arbitrary  frequency  modulation  as  shown 
in  figure  2.  The  pulse  waveform,  p(t),  can  be  expressed  analytically  as: 

(II-l)  j>({)  z  CoS  *■  *?(%>)  *  <Jz] 

for  T  *  t  *■  T _ +TL  where: 

po  po  r 


T  „ 
po 

■  pulse  starting  time 

fc 

■  pulse  carrier  frequency 

PC-Tpo) 

-  initial  pulse  carrier  phase 

kf 

■  constant  of  proportionality 

m(x) 

■  frequency  modulating  function 

-  pulse  width 

The  instantaneous  frequency  function,  f^(t),  associated  with  this  pulse  is: 


(II-2) 


Figure  2  -  Arbitrary  Frequency  Modulation  Function,  m(t) 


The  spectrum  of  the  pulse  waveform  can  then  be  expressed  as  the 


Fourier  Transform  of  p(t)  in  equation  (II-l): 

(n-3)  P(f)  = 


But  since  p(t)  is  nonzero  only  during  the  pulse  duration: 

(II-4)  P(f)  -  ji(t) 


The  difficulty  in  evaluating  this  expression  lies  in  the  fact  that  p(t) 
is  a  function  of  m(t).  Since  m(t),  in  general,  is  not  representable  in 
analytic  form,  its  integral  and  ultimate  application  to  equation  (II-4) 
is  mathematically  untractable.  An  exact  spectrum  for  a  pulse  with  arbi¬ 
trary  frequency  modulation  can  be  generated  using  equations  (II-l)  and 
(11-4)  for  only  certain  m(t).  However,  the  concept  of  approximating  m(t) 
across  the  pulse  holds  promise  for  generating  approximate  spectra. 

Since  integration  is  a  linear  operation,  the  integral  in  equation 
(11-4)  can  be  partitioned  and  rewritten  as  the  sum  of  n  Integrals.  This 
is,  in  essense,  the  decomposition  of  the  waveform,  p(t),  into  n  contiguous 
cells.  The  linear  additivity  property  of  the  Fourier  transform  then  allows 
equation  (11-4)  to  be  written  as: 

ax-5,  m*  ZJ 

<*o  +.  I 

h 

where  in  the  case  at  hand  tQ  ■  Tpo  and  tQ" 

The  spectrum,  then,  of  the  pulse  waveform  with  arbitrary  frequency 
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modulation,  p(t),  can  be  thought  of  as  being  composed  of  the  sum  of  the 
spectra  of  elemental  cells.  If  each  elemental  cell  has  a  characteristic 
spectrum  associated  with  it,  then  P(f),  the  total  pulse  spectrum,  can  be 
generated  quite  easily.  However,  the  problem  of  dealing  with  the  untrac- 
table  form  of  m(t)  across  the  pulse  and  elemental  cells  still  persists. 
There  are  cases,  however,  in  which  the  approximation  of  m(t)  across  an 
elemental  cell  results  in  that  cell  having  a  spectrum  that  is  representable 
in  closed  analytic  form.  For  example,  consider  the  case  where  m(t)  is 
approximated  by  a  constant  value  across  an  elemental  cell.  Let  the  left 
end  point  of  the  cell  be  TQ,  the  cell  width  be  the  initial  phase  of 
the  cell  be  5 t(To)  and  the  radian  center  frequency  of  the  cell  be  w . 

The  characteristic  spectrum  of  this  constant  frequency  unit  cell,  Z(f), 
is  of  the  form  sin(x)/x  and  is  written  as: 


(II-6)  Z(f)  s  ^ 


Thus  the < elemental  cell's  spectrum  can  be  generated  for  any  frequency 
component  as  long  as  the  four  cell  parameters  (TQ,£,  <•>«  and  )  ) 

are  known. 

A  similar  type  of  characteristic  spectrum  results  if  m(t)  is  linearly 
approximated  across  an  elemental  cell.  Consider  the  Instantaneous  fre¬ 
quency  function  shown  in  figure  3.  Let  a  unit  cell  be  selected  as  shown 
extending  from  TQ  to  TQ+  T .  Then  the  Fourier  Transform,  Z(f),  of  this 


elemental  cell  is: 


where  pg(t)  is  Chat  segment  of  p(t)  (equation  II-l)  contained  within 
the  elemental  cell  under  consideration.  The  Instantaneous  frequency  of 
this  cell  is: 

(II-8)  £/t)  =  fc  *  +  t 

Let  the  instantaneous  frequency  (f^Ct))  be  linearly  approximated  (f^(t) ) 
as  shown  in  figure  3,  then: 

m-9,  fa  */(*•*> 


where ^  is  the  frequency  sweep  rate^  *  IT .  Such  linear  FM  is 

called  "chirp"  modulation.  Using  equation  (II-9) ,  the  approximated 
instantaneous  phase  can  be  expressed  as: 

(ii-ioa)  $(t)  5  #&)  y»Cx-T*)]cU 


(II-10b) 


4>(t)  =  4>(To)  +  jft  f,  (t  -  To)  +  1Jp(i -Tit  Z,  ro*r 


This  gives  rise  to  the  approximate  form  of  the  Fourier  Transform  of  the 
elemental  cell. 

(11-11)  2(f)  z Jt 

By  expressing  the  cosine  as  a  sum  of  exponentials,  completing  the  square 
of  the  exponential  arguments  and  letting 
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(II-12a)  2,  s  l/jf  j (i -70)+ 

(II-12b)  £j  =  {?[(*-%)*  Of] 


equation  (II- 11)  reduces  to: 


ui-13)  z(f)--^>\e 


where: 


(II-14a) 


-  'HtJ  -  JirfT0  - 


(II- 14b) 

(II-15a) 

(II- 15b) 

(II-15c) 

(II-15d) 


fi»  5  Ur.) +Jirfl  -jffr.ff 
z,H-.{p[T  *  <£p0J 
hh  i5?\t  *  ^£'j 


The  integrals  in  equation  (11-13)  are  the  Fresnel  integrals.  They  can 
be  evaluated  with  the  use  of  tables,  the  Cornu  spiral  or  the  f'-method 
of  Lanczos(2).  Equation  (11-13)  can  be  written  more  compactly  as: 


(11-16) 


3fr\  .  -X,  C(m)  -C(HL)  +/(s(?I») -S(HL\) ] 

^  CCUL)  ts(uti-s(un1ft\ 
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where  the  Fresnel  Integrals  are  given  by: 


(II-17a)  S(t)  *  J 

JC 

(II-17b)  C(Tf)  =  f  COs(f%  )dv 


Thus  the  Fourier  Transform  of  the  elemental  cell  with  a  linear  frequency 
sweep  can  be  obtained  in  a  straight-forward  manner  by  evaluating  the 
Fresnel  Integrals.  The  derivation  provided  assumes  that  the  frequency 
sweep  within  the  elemental  cell  is  positive.  It  can  be  shown  that  if 
the  sweep  rate  is  negative  the  governing  equation  is  similar  to  equation 
(H-16), 


(11-18) 


Z(f)  -  pp 


C(ZIH)  -C(Ul)  y(sU/if)  - SCul ))] 

C  (tut)  -  C(&tt)-/(S(iJL)’S(VH$]} 


where  the  primed  variables  are  now  defined  as: 
(II-19a)  fi'%  =■  <j>(T0)  ‘J/rfT0  +  *£(/, -rf 

(II-19b)  /ij  -  <#V«) 

(U-20a)  HL '  - 

(II-20b)  2/*'  - 


(II-20c) 


Z1L'  * 


f  i 


(II-20d) 


2JH 


Thus  as  with  the  case  of  a  constant  frequency  elemental  cell,  the  spectrum 
of  a  cell  with  a  linear  frequency  sweep  can  be  generated  for  any  frequency 
component  given  the  cell  parameters  iTtTQt  *To).£i.  and  f  « 

In  summary,  an  approach  has  been  developed  using  an  approximation  to 
arbitrary  FM  to  generate  an  approximation  to  the  pulse  spectrum.  The 
instantaneous  frequency,  f^(t),  is  linearly  approximated  over  an  arbitrary 
number  of  elemental  cells  of  arbitrary  duration.  The  Fourier  transforms 
of  these  elemental  cells  (of  the  constant  frequency  or  chirp  form)  are 
then  coherently  summed  to  produce  an  approximation  to  the  original  pulse 
spectrum.  Appendix  A  is  a  program  listing  for  implementation  of  the 
elemental  cell  model.  Nomenclature  and  notation  in  the  listing  parallels 
that  of  the  development  above.  The  program  in  the  appendix  treats  the 
unit  pulse  problem  as  a  special  case  of  the  more  general  model  developed 
in  section  IV. 


B.  ILLUSTRATION  OF  THE  ELEMENTAL  CELL  MODEL  FOR  THE  UNIT  PULSE. 

It  is  desirable  to  model  an  FM  pulse  for  which  an  exact  spectrum  can 
be  generated  so  that  some  insight  into  and  assessment  of  the  model  per¬ 
formance  can  be  made.  It  is  equally  desirable  to  model  a  pulse  in  which 
the  Instantaneous  frequency  profile  is  dynamic;  l.e.,  constantly  changing 
across  the  pulse.  Both  of  these  objectives  are  achieved  by  selecting  a 
pulse  which  is  subject  to  sinusoidal  frequency  modulation.  Appendix  B  ad¬ 
dresses  the  generation  of  spectra  for  pulses  with  sinusoidal  FM. 
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Consider  Che  case  of  a  5  microsecond  pulse,  which  has  an  instanta¬ 
neous  frequency  function,  f^(t),  given  by: 

(11-21)  J £(t)  COslufrt] 

where 

fc  *  carrier  frequency  *  30  MHz 

kf  *  modulation  index  >1.0 

fm  ■  modulating  frequency  *  200  kHz 

This  instantaneous  frequency  function  is  illustrated  in  figure  4.  The 
figure  demonstrates,  using  two  examples,  the  linearization  and  elemental 
cell  formation  process.  In  case  I,  two  element  cells  are  used  to  generate 
the  spectrum  while  in  Case  II,  six  cells  are  used.  Assume  that  the  desired 
frequency  resolution  is  25  kHz  and  that  the  primary  interest  rests  in 
the  first  three  upper  and  lower  spectral  sidelobes.  In  generating  the 
spectrum  then,  80  frequency  components  are  required  at  25  kHz  intervals 
starting  at  29  MHz.  It  is  worth  noting  that  in  the  elemental  cell  model, 
there  is  no  requirement  that  the  frequency  or  time  intervals  be  of  equal 
size.  Since  the  spectral  components  are  generated  by  evaluating  analytic 
expressions  for  specific  frequencies,  the  spectral  values  for  any  combi¬ 
nation  of  frequencies  can  be  produced. 

Figure  5  displays  the  model-generated  amplitude  spectra  for  the  cases 
developed  in  figure  4.  The  solid  curve  represents  the  exact  amplitude 
spectrum  developed  according  to  Appendix  B.  It  is  not  at  all  unexpected 
that  the  6  cell  model  provides  a  more  accurate  solution  than  the  2  cell 
model.  It  is  of  interest,  however,  that  the  two  cell  model  in  itself 
provides  a  good  approximation  to  the  exact  spectrum.  The  question  that 
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naturally  arises  Is  how  well  do 
those  of  the  actual  spectrum? 


III.  ERROR  ANALYSIS  APPLIED 
TO  THE  ELEMENTAL  CELL  MODEL  FOR  THE 
UNIT  PUI.SE 


A.  ERROR  BOUND 

How  accurate  is  the  elemental  cell  model?  How  many  cells  are 
required  to  achieve  a  given  accuracy?  More  basically,  how  does  the 
model  accuracy  change  as  the  number  of  elemental  cells  is  increased 
or  decreased? 

This  section  assesses  the  elemental  cell  model  performance  in  terms 
of  error  bounds.  The  establishment  of  these  bounds  is  based  on  an 
approach  used  by  Gerald(3)  and  is  developed  in  conjunction  with 
Schoenstadt(4) . 

In  analyzing  the  errors  associated  with  the  elemental  cell  model 
essentially  two  pulses  are  being  compared: 

(1)  the  actual  pulse,  p(t) 

(III-l)  j>(t)  *  COS  L$(tj}  + 

(2)  the  approximated  pulse,  $(t) 


(III-2) 


p(t)  =  COS  [  ^d)]  TpS  t  *  Tp,  +  Tp 


where: 


Tp0  ■  pulse  start  time 

■  pulse  duration 

^>(t)  ■  pulse  instantaneous  phase 

&t)  ”  pulse  approximated  instantaneous 

phase  based  on  the  linearized  instantaneous 
frequency  function 

In  analyzing  the  spectral  error,  the  instantaneous  frequency  error 
is  first  addressed.  This  forms  the  foundation  for  the  phase  error 
analysis  and  ultimately  the  pulse  and  spectral  error  analysis. 

Consider  the  instantaneous  frequency  function,  f^(t),  shown  in 
figure  6.  For  simplicity,  n  elemental  cells  of  equal  length,  h,  can 
be  used  to  generate  the  piecewise  linear  approximation,  £^(t),  such  that: 

(III-3)  f,  (t)  z  >«#-**)] 

where :  fjt)-  linear  approximation  of  the  instantaneous  frequency 

function  across  the  cell. 

/H$x)  m  unit  step  function 

If  it  is  assumed  that  f^(t)  is  at  least  twice  differentiable, 
Gerald (3)  shows  that  the  error  introduced  by  the  linear  approximation  in 
the  cell,  e^(t),  can  be  expressed  as: 

(iu-4)  e.ft)  =  f'lO-fiJt)  . 
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Figure  6  -  Generation  of  f ± (t)  from  f1(t)  Using  Cells  of  Equal  Length 


where  t^j^t^t^  ^  t  k-lA*VA  ck* 

In  the  actual  pulse  equation,  however,  the  Instantaneous  frequency 
error  is  reflected  as  an  error  In  the  Instantaneous  phase;  namely. 


(III-5)  -  kt)  -  ffyo)  '  krpt)  ]js 

Tpo 

Asstime  that  ^CTpo) ,  the  initial  phase,  is  known.  Since  ^(Tpo)  can  be 
set  arbitrarily,  allow  ^(Tpo)  to  be  equal  to  $^(Tpo)  causing  the  first 
two  terms  on  the  right  to  vanish.  Over  the  interval  Tpo  to  Tpo  +  tp  , 
for  any  t  there  exists  a  t^*  TpQ+kh  such  that  t*t^.  Selecting  k  to  be 
the  smallest  integer  such  that  t*tfc,  the  phase  error  then  at  any  time, 
t,  can  be  bounded  from  above  in  the  following  manner: 


(III-6) 


This  integral  can  be  written  as  the  sum  of  integrals  over  the  k  individual 
cells: 


(HI-7) 


*  J"trJ  | 

t:., 


where  -  Tpo  +  jh.  By  construction,  the  limits  of  integration  are  such 
that  the  integrand  in  each  term  of  the  summation  is  given  by  equation 
(III-4) .  Thus, 

tlV 


(III-8) 


for  J-i  v-  By  a  change  of  variable,  let  />W  -  S  - 


so  that 


(HI-9)  ^ /( w(w-i)H^  (n)|jw 


for  /.  The  theorem  of  the  mean  for  Integrals  states  that 


(III-10)  =■ 


provided  that  g(x)  does  not  change  sign  on  the  interval  of  integration 
and  that  \  is  found  in  the  interval  of  integration.  By  virtue  of  this 
theorem,  equation  (III-9)  can  be  written: 

(III-ll)  !#*)-&)( 14^1  JiwCw-oldw 

for  OSiy&l  .  Equation  (III- 11)  essentially  expresses  the  bound  on  the 

instantaneous  phase  error  as  the  accumulation  of  errors  over  k  individual 

\ 

elemental  cells.  An  upper  bound  can  be  set  on  the  second  derivative  of 
the  instantaneous  frequency  in  each  cell  such  that  for  the  cell  the 

bound  is  Mj  -  max  for  tjm%  .  Making  use  of  this  bound  and 

Integrating  equation  (III- 11)  yields: 

(III-12)  J^t)  -  #*)|  “  ^ 

More  generally  the  second  derivative  can  be  bounded  over  the  entire  pulse 
such  that  M  ■  max  Mj  for  j  ■  1,2,3 . ,k.  It  follows  then  that 


equation  (III-12)  becomes: 


(III-13) 


j  *>-*>/*  ^ 


The  original  selection  of  k,  however,  required  that  T  +h (k-1 ) 4 1 iT  +hk . 

po 

In  other  words,  kh*t-Tp0+h.  Making  use  of  this  and  the  fact  that  hn  mfp, 
equation  (111-13)  can  be  written: 

(III- 14)  ±  [t-Tp  *5] 


Thus  an  upper  bound  on  the  absolute  phase  error  has  been  obtained  as  a 
function  of  time.  It  now  remains  to  relate  this  error  to  the  actual 
pulse  error. 

The  absolute  pulse  error  can  be  formed  using  equations  (III-l)  and 
(III-2)  as: 


(III-15) 


(III-16) 


=  |  C  03  [<£(0]  -  cosDftt)]  | 


Since  the  sin(x)  is  less  than  or  equal  to  one, 
(III-17)  |  p(t)  -  p(t)  |  ±  jjS/AJ 
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By  making  use  of  the  fact  that  lsln(x)J  *  |x| ,  the  pulse  error  bound  reduces 


to: 

(111-18) 

dii-19)  |  pft)  -  pit)  |  ±  £0e  (j .  ^ ^ 


With  this  knowledge  of  the  pulse  error,  the  bound  on  the  spectral  error 
can  be  addressed.  Writing  the  spectral  error  In  terms  of  the  pulse 
error, 

(III- 20)  p(f)  -P(f)  s  J  [ 

Tpo 


where  P(f)  and  f’(f)  are  the  exact  and  approximated  spectra  respectively. 
The  absolute  error  can  be  written: 


(III-21) 


\pcf)-P(f)\i  jTp«)-p«U  Jt 

Tpm 


Using  the  derived  expression  for  the  absolute  phase  error  (equation 


III-14) , 

(in-22)  \W)-fo)\  i  ' Tr° 

TP» 

(HI-23)  |p(f)-fe)|  * 


Additionally,  since  the  absolute  value  of  a  difference  Is  greater  than  or 
equal  to  the  absolute  value  of  the  difference  of  the  absolute  values,  the 


bound  on  the  spectrum  error  using  the  cell  model  can  be  expressed  as: 

«n-*>  jjm/-lkf)lj  -  00 f  i ) 

In  general  then,  the  accuracy  of  the  elemental  cell  model  depends  on 
the  ma-g-tumm  value  of  the  second  derivative  of  the  instantaneous 
frequency  function  (a  measure  of  how  nonlinear  the  modulation  is) , 
the  pulse  width  and  the  number  of  elemental  cells  which  subdivide  the 
pulse.  It  is  significant  to  note  that  the  model  converges  to  the  actual 
spectrum  inversely  as  the  square  of  the  number  of  cells. 

While  the  error  bound  developed  above  assumed  elemental  cells  of 
equal  length  and  bounded  the  second  derivative  of  the  Instantaneous 
frequency  function  across  the  entire  pulse,  a  tighter  bound,  and  one 
which  is  more  generally  applicable  to  elemental  cells  of  unequal  length, 
can  be  obtained  by  bounding  the  error  contribution  from  each  Individual 
cell  and  then  summing  these  contributions  over  the  total  number  of 
elemental  cells.  The  total  spectrum  error  is  given  by 

(in-25)  pen  * p(f )  *  Efi  w  + + . 

where  P1(f)  and  P^(f )  are  the  true  and  approximated  spectra  respectively 
of  the  i*-*1  elemental  cell.  So  that  the  absolute  spectrum  error  can  be 
bounded  from  above  by  summing  over  the  absolute  error  associated  with 
each  elemental  cell. 

(HI-26)  / P(f)  - P(f)\  *  £[  |  ? (f)  “  ft (0 1 

1  fl 
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The  absolute  spectrum  error  associated  with  each  elemental  cell, 
lP(f)-%f)l  ,  remains  to  be  developed.  Consider  the  kc^  elemental 
cell  as  depicted  in  figure  7.  An  expression  analogous  to  equation 
(III-5)  can  be  developed  to  describe  the  absolute  instantaneous  phase 
error  at  some  time,  t,  associated  with  the  kc^  elemental  cell. 


(III— 27) 


jm  -  kt)j  i  jrt'o  -  &u|  *  j*J  |  £.6) 


where  &  t  -  Toh*^k  —  7a(/[+i).  In  addressing  the  analogous  equation 
(III-5),  the  phase  error  at  the  beginning  of  the  pulse  was  set  equal 
to  zero  since  ^(T0)  was  arbitrarily  set  equal  to  ^(T0) .  In  the  general 
case,  in  which  successive  elemental  cells  are  being  considered,  &T„k> 
is  not  arbitrary.  The  phase  difference  at  Tq^  must  be  viewed  as  a 
cummuiative  error  resulting  from  the  approximations  made  in  the  succes¬ 
sive  elemental  cells  occurring  before  T^.  To  develop  this  point,  assume 
as  before  that  at  the  beginning  of  the  pulse  (that  is  the  first  elemental 
cell)  that  Using  the  above  equations,  the 

phase  error  at  the  end  of  the  first  cell  or,  equivalently,  at  the  begin¬ 
ning  of  the  second  cell  is: 


(III-28) 


16m- 


Using  the  expression  for  the  elemental  cell  instantaneous  frequency  error 
developed  in  equation  (III-4) , 

Jo,  t  f, 

(III-29) 


for  Tol*Hj*T0|  +t( .  Applying  the  mean  value  theorem  of  integrals, 
bounding  the  second  derivative  of  the  instantaneous  frequency  across 
the  cell  by  its  maximum  value,  M^,  and  integrating  yields: 

(III-30) 


Thus  the  phase  error  at  the  beginning  of  the  second  cell  is  a  function 
of  the  quality  of  the  approximation  made  in  the  first  cell.  A  similar 
development  can  be  applied  to  generate  the  phase  error  at  the  end  of 
the  second  cell  (beginning  of  the  third  cell) ;  the  absolute  phase  error 
at  the  beginning  of  the  third  cell  being: 


(111-31) 


The  first  term  on  the  right  is  given  by  equation  (HI- 30)  and  the  integral 
can  be  evaluated  as  above  to  yield: 


(III-32) 


Extending  this  result,  the  absolute  phase  error  at  the  beginning  of  the 
kth  cell  can  be  expressed  as: 

x-i  j 

(HI-33)  /#£*)- &5a)|  £  £ 

This  result  can  be  applied  to  equation  (III-27)  along  with  the  expression 
for  the  instantaneous  frequency  error  in  equation  (III-4)  to  develop  an 
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expression  for  the  absolute  phase  error  associated  with  the  kch  cell 
as  a  function  of  time: 


(III- 34) 


for  T0k*ttT0k+?<.  This  expression  for  the  phase  error  can  be  used  in 
combination  with  equation  (III-21)  to  form  the  upper  bound  on  the 
spectral  error  associated  with  the  kch  elemental  cell. 

(III- 35)  J PK(f)-?Jf)l  -  lj 


(III-36) 


This  represents  then  a  bound  on  the  local  spectrum  error  associated  with 
the  kch  elemental  cell.  To  obtain  a  bound  on  the  global  error  across 
the  pulse,  the  local  error  of  each  cell  can  be  summed  over  all  of  the 
elemental  cells,  as  in  equation  (III-26),  so  that 


(III-37) 


liPrt-Mji  £  Or  1^9 


,  tfAfX 

/ JL 


While  this  expression  for  the  error  bound  is  more  unwieldy  than  that 
developed  in  equation  (III-24) ,  it  has  provisions  for  addressing  ele¬ 
mental  cells  of  non-uniform  length.  It  also  bounds  the  second  derivative 
of  the  instantaneous  frequency  function  on  a  cell  by  cell  basis,  yielding 
a  tighter  spectral  error  bound.  This  tighter  bound  is,  of  course,  at  the 


expense  of  complexity  of  evaluation.  The  determination  of  the  bound 
on  the  second  derivative  of  the  instantaneous  frequency  function,  M, 
for  equation  (111-24)  may  be  difficult  enough.  To  further  determine 
the  bounds  Mj  for  all  elemental  cells  for  application  to  equation 
(111-37)  may  be  unacceptably  tedious.  A  trade-off  can  be  made  to  loosen 
the  error  bound  of  equation  (111-37)  while  preserving  the  provision  for 
addressing  non-uniform  elemental  cell  lengths  by  bounding  the  second 
derivative  of  the  instantaneous  frequency  function  across  the  pulse 
(i.e.,  M  »  max  ,  j  -  1,2, . n) 

an- 38)  J-M 1  *  *#•) 

Thus  various  approaches  of  different  degrees  of  compexity  can  be 
used  to  develop  theoretical  maximum  error  bounds  on  the  spectrum  error 
associated  with  the  elemental  cell  model.  The  choice  of  approach  must 
be  based  on  problem  complexity  and  a  priori  knowledge  of  the  pulse 
structure.  The  most  important  insight,  however,  that  the  maximum  error 
bound  provides  is  that  in  applying  the  elemental  cell  model,  this  error 
tends  toward  zero  inversely  as  the  square  of  the  number  (n)  of  cells. 

B.  EMPIRICAL  ASSESSMENT 

The  objective  of  this  section  is  to  assess  the  actual  elemental 
cell  model  performance  in  the  light  of  the  error  bounds  developed  in 
section  III-A.  To  achieve  this  a  unit  pulse  with  sinusoidal  FM  is 
examined  (as  in  Section  II-B) .  The  pulse  is  such  that  the  carrier 


frequency,  fc,  is  30  MHz;  Che  modulating  frequency,  fm,  is  200  kHz; 
che  modulation  index,  kf,  is  one;  and  che  pulse  duration,^,  is  five 
microseconds.  The  exact  spectrum  of  this  pulse  can  be  determined  (see 
Appendix  B) .  Using  the  elemental  cell  model  (Section  II-B) ,  approximate 
spectra  can  be  generated  for  cases  where  the  number  of  elemental  cells, 

n,  is  1,2,3, .  Each  of  these  approximated  spectra  can  then  be 

compared  to  the  exact  spectrum  to  determine  the  maximum  spectrum  error 

calculated  for  n=l,2,3 .  (For  the  purpose  of  this  analysis, 

frequency  components  between  29  MHz  and  31  MHz  are  examined  at  20  kHz 
intervals . ) 

Since  the  unit  pulse  under  consideration  has  an  instantaneous 
frequency  function,  f^(t),  given  by 

(III-39)  ~  t  +  ffft*  C°S 

the  maximum  bound  on  the  second  derivative  of  the  instantaneous  fre¬ 
quency  function  can  be  determined  analytically.  Since  this  is  the  "M" 
of  equation  (III-24),  that  equation  can  be  evaluated  to  yield  the  maximum 
theoretical  spectrum  error  bound  as  a  function  of  n,  the  number  of 
elemental  cells  used  to  approximate  the  pulse. 

It  is  also  possible  to  generate  a  tighter  theoretical  bound  on  the 
spectrum  error  by  using  equation  (III-37)  and  determining  from  equation 
(III-39)  the  bound  on  the  second  derivative  of  the  Instantaneous 
frequency  function  across  each  Individual  cell. 

Thus  a  set  of  three  curves  can  be  generated: 

(1)  The  maximum  theoretical  spectrum  error  bound  with 


bounded  across  the  pulse  (equation  III- 24)  as 
a  function  of  n,  the  number  of  elemental  cells; 

(2)  The  maximum  theoretical  spectrum  error  bound  with  f*"(t) 
bounded  across  each  elemental  cell  (equation  III-37) 

as  a  function  of  n; 

(3)  The  maximum  calculated  spectrum  error  as  a  function  of  n. 

Figure  8  presents  this  set  of  curves  for  the  pulse  with  sinusoidal  FM 
discussed  above.  The  first  theoretical  bound  (equation  III-24)  is 
represented  by  the  solid  line;  the  second  (equation  III-37),  by  the 
dashed  line;  and  the  calculated  error,  by  the  dotted  line.  The  tightening 
of  the  maximum  theoretical  error  by  bounding  the  second  derivative  of 
the  instantaneous  frequency  function  across  each  elemental  cell  vice 
across  the  pulse  can  be  observed  quantitatively  by  comparing  the  solid 
and  dashed  curves.  The  maximum  calculated  error  (dotted  curve)  is  well 
within  the  bounds  set  by  equations  (III-24)  and  (III-37).  In  fact,  the 
calculated  error  is  generally  an  order  of  magnitude  smaller  than  the 
maximum  error  bound  set  by  equation  (III-24).  Inspection  of  Figure  8 
indicates  that  all  three  curves  generally  follow  the  characteristic 
convergence,  which  goes  inversely  as  the  square  of  the  number  (n)  of 
elemental  cells. 

An  important  consideration  in  applying  the  elemental  cell  model  is 
the  determination  of  the  number  of  elemental  cells  to  be  used.  Obviously, 
the  greater  the  number  of  cells,  the  closer  the  approximated  spectrum  is 
to  the  exact  spectrum,  but  at  the  cost  of  complexity  and  processing  time. 
The  establishment  of  the  inverse  squared  relationship  between  the  number 
of  cells  and  the  theoretical  maximum  error  bound  allows  the  determination 


Figure  8  -  Theoretical  Error  Boiind  and  Maximum  Calculated 
Error  as  a  Function  of  the  Number  of  Elemental 


Cells  Used  in  the  Model  Process 


of  the  minimum  number  of  cells  required  as  a  function  of  the  maximum 
tolerable  error.  Let  €  be  the  maximum  tolerable  error,  then  equation 
(III-24)  can  be  written  as: 


(111-40) 


so  that 

(III-41) 


or  for  the  case  when  the  number  of  cells  in  large: 


(III-42) 


For  example,  in  the  case  of  the  unit  pulse  with  a  200  kHz  sinusoidal 
FM  discussed  above,  the  spectrum  has  a  peak  value  on  the  order  of 
(^»/2)  2.5x10“6.  if  the  interest  is  primarily  in  the  mainlobe,  the 
maximum  tolerable  error  may  be  an  order  of  magnitude  less  than  the  peak 
value,  so  that 2.5x10”^.  Then,  knowing  the  value  of  M  and  £, 
equation  (III-42)  can  be  evaluated  (or,  equivalently  Figure  8  could  be 
consulted)  to  determine  that  15  cells  (n  -  14.3)  are  required  to  assure 
that  the  maximum  tolerable  error  is  not  exceeded.  In  this  example, 
n»15  provides  considerably  more  accuracy  than  is  required,  but  does 
also  provide  the  guarantee  that  the  error  will  not  be  greater  than  the 
accepted  limit. 

In  practice,  the  second  derivative,  used  freely  above,  may  seldom  be 
available  in  analytic  form.  In  those  cases,  a  sizing  of  the  problem  can 


w 


be  made  using  Che  central  difference  approximation  to  the  second  deriv¬ 
ative  (3).  Specifically,  given  a  function  g(t),  as  in  figure  9,  the 
approximation  to  the  second  derivative  at  TQ  is: 


(III-43) 


In  the  case  of  the  test  signal  discussed  above,  the  maximum  of  the 
second  derivative  of  the  instantaneous  frequency  function  across  the 
pulse  is  given  by:  M  ■  3.1583x10*7.  if  the  central  difference 
approximation  of  equation  (III-43)  is  applied  to  the  test  pulse  at 
T0  -  2.5  microseconds  for  various  values  of  h  and  the  resulting 
approximated  derivatives  applied  to  equation  (III-42) ,  the  following 
estimates  of  the  numbers  of  cells  required  to  produce  an  error  no 
worse  than  2.5xlO-7  are  obtained: 


for  h  ■  1.2500/ts;  f^(T0)  *  2.56x10*7;  n  *  12.6 
for  h  -  0.6250 >3}  ?"(T0)  *  3.00xl017;  n  -  13.7 
for  h  -  0.3125 -  3.12xl017;  n  -  13.9 


Thus  this  approximate  approach  to  sizing  the  modeling  problem  is  in 
close  agreement  with  the  exact  solution  of  n  *  14.3. 

The  maximum  error  bounds  developed  in  Section  (III-A)  provide  a 
valuable  insight  into  the  nature  of  the  convergence  of  the  elemental 
cell  model  and  a  useful  tool  in  determining  the  number  of  elemental 
cells  required  to  model  a  specific  pulse. 
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C.  ASSESSMENT  OF  ERRORS  ACROSS  THE  SPECTRUM 


To  this  point,  the  error  analysis  has  addressed  the  worst  case, 
maximum  error  bound,  aspects  of  the  elemental  cell  model  and  has 
emphasized  the  1/n^  convergence.  But  this  maximum  error  is  simply  a 
bound  and  does  not  address  the  error  associated  with  individual  fre¬ 
quency  components  across  the  spectrum. 

Consider  the  difference,  E*  ,  between  the  exact  and  the  approximated 
pulses  in  the  time  domain, 

(III-44)  (t)  =  Co  sCAfcfl  -  CosC^O] 

where  $t)  is  derived  from  the  instantaneous  frequency  function,  £j(t), 
equation  (II-2)  and  At)  is  similarly  derived  from  the  linearly  approx¬ 
imated  instantaneous  frequency  function,  f^(t).  Then  At)  and  At)  are 
related  by: 

(III-45)  <P(t)  *  4>(t)  *  6(t) 


where  6(t)  is  the  phase  error  as  a  function  of  time  such  that 

(III-46)  e(i)  <PH)-Wt)  -  Jnjtfan-fatx Jdx 

e 

(III-47)  6(t)  =  /eu)<4* 


where  e(t)  is  the  instantaneous  frequency  error  as  a  function  of  time. 


Substituting  equation  (III-45)  into  (III-44)  yields: 


(111-48)  £t  (i )  -  COS  L$(t)  *  -  COS  Iftti 


(III-49) 


£i({)  -  JLs/fj  ($it)  * 


If  tf(t)  is  small,  that  is,  for  example,  if  £(t)  <•  0.5  (arbitrarily), 

then 


(III-50)  ^  ^  *  <ft)  Slid  f$(i)  + 

Additionally,  if  then 

(in-51)  £t(i)  *  e(i)s,A/l<ht )] 

If  <(t)  has  a  Fourier  Transform,  £(f),  and  sin^(t)  has  a  Fourier  Trans¬ 
form,  £(f);  then  the  error  associated  with  Et(t)  in  the  frequency  domain, 
Ef(f),  is  given  by  the  convolution  of  £ (f)  and  1*(f ) . 

(III-52) 

While  this  egression  is  far  from  simple,  it  may  provide  an  insight  into 
the  behavior  of  errors  across  the  spectrum  as  the  number  of  elemental 
cells  is  varied. 

Consider,  again,  the  teat  pulse  of  unit  amplitude  with  sinusoidal  FM. 
If  this  pulse  is  modeled,  using  four  elemental  cells,  the  instantaneous 
frequency  error,  e(t),  and  the  instantaneous  phase  error,  tf(t),  can  be 


generated  as  in  figure  10.  Before  applying  the  instantaneous  phase 
error  function  to  equation  (III-52),  let  us  first  examine  the  appli¬ 
cability  of  the  equation  to  the  case  at  hand. 


In  arriving  at  equation  (111-51)  and,  ultimately,  (III-52),  two 
assumptions  were  made:  (1)  that  the  instantaneous  phase  error, €(t), 
was  small  (that  is  that  tf(t)a  sin(e(t))  and  (2)  that  $(£)  • 

Figure  10  illustrates  that  assumption  one  is  satisfied;  i.e.  6(t)*-0.5. 
As  for  assumption  (2),  generally,  J^(t)  behaves  as  &t),  so  that 


(III-53) 


x&OfrT-  S/aj  (o.VttT) 


for  T  in  microseconds.  Realizing  that  this  function  is  based  on  modulo 
two  pi,  there  will  be  a  few  cases  where  <(t) ~  ,  but  in  general  over 

jk  * 

the  pulse  ^Ct)»  .  (See  Figure  11). 

Consider  now  the  instantaneous  phase  error,  £(t),  in  figure  10,  in 
the  light  of  equation  (III-52).  The  Fourier  transform  of  €(t)  is 
essentially  a  sin(x)/x  function  centered  about  the  carrier  frequency, 
in  this  case  200  kHz.  In  actuality,  the  spectrum  is  not  quite  this 
straightforward,  because  of  the  relative  proximity  of  carrier  frequency 
to  the  origin  (f«0) .  The  spectrum,  as  shown  in  figure  12,  is,  in  general 
characterized  by  two  main  lobes  of  400  kHz,  sidelobes  of  200  kHz,  and 
a  spacing  between  the  mainlobes  on  the  order  of  440  kHz. 

The  basic  form  of  the  spectrum,  't(f),  generated  by  the  instantaneous 
frequency  function,  ^(t) ,  is  of  the  form  shown  in  figure  13.  Ef(f), 
can  be  considered  generally  based  on  figures  12  and  13.  It  is  expected 
that  Ef(f)  would  have  a  basic  sin(x)/x  form  with  two  mainlobes  being 
about  400  kHz  wide  and  symmetrically  spaced  220  kHz  above  and  below  the 


Figure  10  -  Time  Domain  Error  Functions  for  a  Four 


Cell  Model  of  a  Pulse  with  Sinusoidal  FM 


Figure  11  -  Comparison  of  Phase  and  Phase  Error  across 
a  Sinusoidally  Modulated  Unit  Pulse  Using 
a  Four  Elemental  Cell  Model;  Percent  of  Time 
that  the  Ratio,  R,  is  Greater  Than  a  Given 
Threshold  Value,  x,  as  a  Function  of  the 


Threshold  Value 
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Frequency [MHz] 


Figure  12  -  Amplitude  Spectrum  of  a  5.0£s  Pulse  with  a 


200  kHz  Carrier 


Spectrum  , 
Amplitude (xlO  ) 


Figure  13  -  Exact  Amplitude  Spectrum  of  a  5^s  Pulse  with 
200  kHz  Sinusoidal  FM  (Modulation  Index  ■  1.0) 
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peak  frequency  component  of  $(f).  Further,  it  would  be  expected  that  a 
negligible  error  would  be  associated  with  the  peak  value  of  £(f)  and 
that  the  upper  and  lower  sidelobes  of  Ef(f)  would  be  about  200  kHz 
wide.  Figure  14  is  helpful  in  addressing  this  qualitative  analysis  of 
Ef(f).  The  solid  line  in  the  figure  represents  the  exact  spectrum  of 
a  unit  pulse  with  a  200  kHz  sinusoidal  frequency  modulation  imposed 
upon  it.  (i.e.  The  test  signal  discussed  above).  The  additional 
projections  (x’s,  circles  and  triangles)  represent  the  spectrum  error, 
E^(f),  associated  with  the  modeling  of  the  test  signal  using  four, 
eight  and  sixteen  elemental  cells,  respectively.  For  the  moment,  the 
case  of  four  elemental  cells  will  be  addressed.  Note  that  the  mainlobes 
of  Ef(f)  are  on  the  order  of  400  kHz  wide  and  that  these  lobes  are 
centered  fairly  symmetrically  about  the  peak  value  of  P(f),  being  offset 
by  200  to  250  kHz.  As  predicted,  the  upper  sidelobes  tend  to  decrease 
in  width  to  the  order  of  200  kHz.  Thus,  while  equation  (III-52)  is 
fairly  unwieldy,  a  basic  insight  into  the  behavior  of  Ef(f)  can  be 
obtained. 

Consider  now  the  two  additional  cases  presented  in  figure  14  for 
the  cases  when  eight  and  sixteen  cells  are  used.  The  Instantaneous 
frequency  and  phase  error  for  an  eight  cell  and  a  sixteen  cell  simulation 
are  shown  as  functions  of  time  in  figures  15  and  16,  respectively.  It 
can  be  seen  that  in  both  cases  the  instantaneous  phase  error  departs 
slightly  from  that  examined  under  the  case  n*4,  but  that  the  general  form 
of  £(t)  persists.  In  the  frequency  domain,  this  is  reflected  in  the 
similar  forms  of  Ef(f)  for  n»4,8,16.  Of  particular  interest  is  the 
scaling  associated  with  each  of  the  three  phase  error  functions.  All 
three  instantaneous  phase  error  functions  are  forms  of  a  perturbed 
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(b)  Instantaneous  Phase  Error  vs.  Time 


Figure  15  -  Time  Domain  Error  Functions  for  an  Eight 
Cell  Model  of  a  Pulse  with  Sinusoidal  FM 
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(a)  Instantaneous  Frequency  Error  vs.  Time 


(b)  Instantaneous  Phase  Error  vs  Time 


Figure  16  -  Time  Domain  Error  Functions  for  a  Sixteen 


Cell  Model  of  a  Pulse  with  Sinusoidal  FM 


200  kHz  waveform  with  varying  amplitudes: 

for  four  cells,  the  amplitude  is  .214  radians 
for  eight  cells,  the  amplitude  is  .051  radians 
for  sixteen  cells,  the  amplitude  is  .012  radians. 

So  that  in  increasing  the  number  of  cells  from  four  to  eight,  the  peak 
phase  error  is  reduced  by  a  factor  of  0.24.  An  increase  in  the  number 
of  cells  from  eight  to  sixteen  reduces  the  peak  phase  error  by  a  factor 
of  0.24.  Considering  then  the  convolution  of  equation  (III-52)  and 
assuming  that  P(f)  is  essentially  constant  for  n-4,8,16,  it  would  be 
expected  that  Ef(f)  would  be  improved  by  about  a  factor  of  0.24  for 
increases  in  the  number  of  elemental  cells  from  four  to  eight  and  from 
eight  to  sixteen.  This  is,  in  fact,  the  case  generally  portrayed  in 
figure  14.  Typically,  the  improvement  experienced  in  going  from  four 
to  eight  cells  and  from  eight  to  sixteen  cells  is  on  the  order  of  a 
0.20  to  0.25  reduction  in  the  spectrum  error.  Again,  this  reduction 
in  spectrum  error  by  a  factor  of  4  when  the  number  of  elemental  cells 
are  doubled  is  characteristic  of  the  1/n^  convergence  discussed  in 
section  III-A.  While  some  abberatlons  occur  in  the  convergence  to  the 
exact  spectrum  (e.g.  between  29.0  and  29.4  MHz),  it  is  significant  that 
the  convergence  across  the  spectrum  is  generally  uniform.  Even  in  cases 
of  anomalies  (29.0  to  29.4  MHz),  the  convergence  tends  to  be  uniform 
as  n  increases. 

In  summary,  by  examining  the  phase  error,  it  is  possible  to  project 
the  basic  form  of  the  spectrum  error,  Ef(f)  and,  perhaps  of  more  signif¬ 
icance,  to  examine  the  behavior  of  the  spectrum  error  as  the  number  of 


cells  Is  varied.  This  provides  an  Important  Insight  Into  E^(f )  and 
the  convergence  of  ^(f)  to  P(f). 


of  the  pulse  s(t)  in  equation  (IV-1)  is  given  by  the  convolution  of 
V(f)  and  the  Fourier  transform  of  cos  ^(t) : 


(IV-3)  s(f)  -  m  «  fit*  =  W) «  Rf) 


where  P(f)  is  given  by  equation  (II-3).  Using  the  approach  of  equation 
(II-5),  the  spectrum,  S(f),  can  be  represented  by  the  sum  of  the  spectra 
of  elemental  cells: 


where  t0*Tp0  and  tn  ■  Tpo+£.  The  spectrum,  2(f).  of  any  individual 
element  cell  is  given  by 


where  Tq  is  the  start  time  of  the  elemental  cell  and  t  is  the  pulse 
duration  of  the  elemental  cell.  Section  II  treated  the  special  case 
of  equation  (IV-5)  when  v(t)  is  equal  to  one  across  the  elemental  cell 
and  developed  equations  (II-6  and  11-13)  for  the  characteristic  elemental 
cell  spectra  for  constant  and  chirp  frequency  cells.  If  the  pulse 
envelope,  v(t),  is  not  equal  to  one,  but  rather  some  constant  across 
the  elemental  cell,  then  equations  (II-6  and  11-13)  still  apply  and 
need  only  be  modified  by  a  scale  factor.  If,  on  the  other  hand,  v(t) 
is  not  constant  across  the  cell,  equations  (II-6  and  11-13)  must  be 
modified  to  account  for  this  non-constant  behavior. 


Consider  the  instantaneous  frequency  function,  f^(t) ,  and  amplitude 


modulation,  v(t),  depicted  in  Figure  17.  In  the  figure,  an  elemental 
cell  extending  from  T0  to  Tq  +  T  is  generated  by  linearly  approximating 
f^(t)  and  v(t)  across  the  elemental  cell.  The  linearized  instantaneous 
frequency  function,  £^(t),  has  been  discussed  in  great  detail  in  Section 
II.  The  linearized  amplitude  modulation,  v(t),  is  given  by: 


(IV-6)  -  A+hU'V 


r 


where  "A"  is  the  initial  cell  amplitude  and  "h"  is  the  slope  of  the 
linearized  envelope.  This  allows  the  element  cell  spectrum, 2(f) »  of 
equation  (IV-5)  to  be  approximated  by: 

JA 

a  I  r -S&jt/j.  .  Li 

(IV-7) 


j9*t  j+t 

Z(f)  =  A +hj(t 

To  To 


The  first  integral  is  simply  the  problem  addressed  in  Section  II 
modified  by  a  multiplicative  constant.  The  second  integral  adds  a 
degree  of  generality  to  the  elemental  cell  model  by  accounting  for  the 
fact  that  the  amplitude  modulation  is  not  constant.  In  consonance  with 
the  model  development  of  Section  II,  there  are  two  special  cases  of 
equation  (IV-7)  that  merit  attention: 


(1)  h^o  and  the  instantaneous  frequency  function  is  constant 
across  the  cell;  and 

(2)  hj*o  and.  the  instantaneous  frequency  function  is  linear 
across  the  cell. 


In  the  first  case,  the  approximated  instantaneous  phase  function,  ^(t) 


is: 


(lv-8)  4>(t)  -  +J#fct 

for  •  By  applying  this  function  to  equation  (IV-7), 

straightforward  integration  yields  the  following  expression  for  the 
characteristic  spectrum  of  an  elemental  cell  with  linear  amplitude 
modulation  and  a  constant  carrier  frequency: 


In  the  generalized  elemental  cell  model,  this  expression  is  the  analog 
of  equation  (II-6) . 

In  addressing  the  special  case  of  equation  (IV-7)  in  which  both 
the  elemental  cell  amplitude  and  frequency  are  linearly  swept,  an 
expression  analogous  to  equation  (11-13)  will  be  developed.  In  this 
case,  the  approximated  instantaneous  phase  function,  p(t),  is  given  by 
equation  (II-10b) .  The  first  Integral  of  equation  (IV-7)  is  simply  a 
constant.  A,  times  an  integral  evaluated  in  Section  (II-A);  namely, 
equation  (11-13).  Using  Euler's  equations,  the  second  Integral  of 
equation  (IV-7)  can  be  arranged  as  the  sum  of  two  integrals;  the  first 
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is  a  function  of  (f^-f )  while  the  second  is  a  function  of  (f^+f ) . 
Designating  the  first  of  these  integrals  by  I^(f)  and  the  second  by 
^(f),  I^(f)  can  be  written: 

7^T 

(iv-io)  !,((■)■-  J(t-r.)eir  jt 


By  completing  the  square  of  the  integrand's  exponential  argument,  and 
defining  z^  as  in  equation  (11-12) ; ,  as  in  equation  (11-14) ;  and 
Z1L  and  Z1H  as  in  equation  (11-15) ,  I-,  (f)  can  be  written: 


The  first  integral  is  directly  integrable.  The  second  is  of  the  Fresnel 
form.  Defining  the  Fresnel  integrals,  C(x)  and  S(x),  as  in  equation 
(11-17),  I^(f)  can  be  written  as: 


(IV-12) 


«»•  $£{<*** -s** 

-  -am)y @p")  -  stvM 


A  simillar  analysis  can  be  applied  to  the  integral  12(f). 

For  simplicity,  define  the  evaluated  Fresnel  integrals  as: 


(iv-i3a)  Ft  *  CCm)  -C(hl)  -  SChl)] 


(IV-13b) 


Fa  =  C(2JH)~C(2JL)  tj[s(ijL)  -5{ZJX)] 


With  this  definition  and  letting  the  limits  of  integration  be  defined 
as  in  equation  (11-15)  and  letting^  and  be  defined  as  in  equation 
(11-14),  equation  (IV-7)  can  be  written  in  integrated  form  as: 


This  expression  is  the  analog  of  equation  (11-13)  for  the  case  where 
arbitrary  amplitude  modulation  is  considered.  This  is  apparent  in  the 
terms  that  have  been  added  that  are  a  function  of  the  envelope's  slope. 

Thus  any  pulse  with  an  arbitrary  amplitude  and  frequency  modulation 
can  be  modeled  by  expressing  the  pulse  envelope  and  the  instantaneous 
frequency  as  linear  piecewise  continuous  functions.  These  linearized 
function  are  then  segmented  into  a  series  of  successive  elemental  cells 
such  that  each  cell  can  be  classified  as  one  of  the  following: 

(1)  linear  frequency  sweep;  linear  envelope 

(2)  linear  frequency  sweep;  constant  envelope 

(3)  constant  frequency  sweep;  linear  envelope 

(4)  constant  frequency  sweep;  constant  envelope 

Associated  with  each  of  the  four  classes  of  elemental  cells  is  a 
characteristic  spectrum,  so  that  once  the  elemental  cell's  class  has 
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been  determined,  Its  spectrum  can  be  generated  directly  from  cell 
parameters.  Finally,  the  total  pulse  spectrum  Is  obtained  by  summing 
over  the  spectra  of  all  elemental  cells. 

The  generalized  elemental  cell  model  Is  flexible  and  easy  to  use. 
There  is  no  requirement  that  the  elemental  cells  be  of  uniform  length. 
In  fact,  judicious  construction  of  the  element  cells  allows  for  the 
possibility  of  dynamic  sampling.  That  is  smaller  elemental  cells 
(higher  sampling  rates)  may  be  used  during  transient  portions  of 
the  pulse  while  larger  cells  (lower  sampling  rates)  may  be  used  during 
more  static  portions  of  the  pulse.  Additionally  frequency  resolution 
is  not  a  function  of  the  sampling  rate  and  the  number  of  samples  used. 
Arbitrary  frequency  resolution  can  be  achieved  since  the  model  uses 
an  analytic  approach  to  spectrum  generation.  For  the  same  reason,  an 
arbitrary  number  of  spectrum  frequency  component  can  be  generated. 

A  program  listing  which  implements  the  generalized  elemental  cell 
model  is  provided  in  Appendix  A. 


B.  ILLUSTRATION  OF  THE  GENERALIZED  ELEMENTAL  CELL  MODEL 


To  illustrate  the  application  of  the  generalized  model,  consider  as 
in  Section  (II-B)  the  sinusoidally  frequency  modulated  pulse  but  allow 
the  pulse  envelope  to  also  be  sinusoidally  shaped  as  shown  in  Figure  18. 
The  spectrum  of  this  pulse  is  discussed  in  Appendix  B.  Let  the  pulse  be 
such  that  the  carrier  frequency,  fc  is  30  MHz;  the  frequency  modulating 
frequency,  fm,  is  200  kHz;  the  frequency  modulation  index,  kf,  is  one; 
the  amplitude  modulating  frequency  is  100  kHz;  and  the  pulse  duration 
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is  five  microseconds.  (Assume,  again,  that  a  25  kHz  resolution  is 
desired  and  that  the  primary  Interest  is  in  the  first  three  upper 
and  lower  spectral  sidelobes.) 


Figure  18  demonstrates  the  generation  of  elemental  cells  for  this 
pulse.  Two  cases  are  considered.  In  case  I,  two  cells  are  used  to 
generate  the  spectrum;  while,  in  case  II,  eight  cells  are  used.  Figure 
19  displays  the  amplitude  spectra  generated  by  these  two  cases.  The 
solid  line  represents  the  exact  spectrum,  the  development  of  which  is 
addressed  in  Appendix  B.  The  general  results  are  not  unlike  those 
obtained  for  the  unit  pulse  in  Section  (II-B) .  In  the  following  section 
it  will  be  shown  that  the  convergence  of  the  generalized  model  parallels 
that  of  the  unit  pulse  model  being  inversely  related  to  the  square  of  the 
number  of  elemental  cells  used  in  generating  the  spectrum. 


V.  ERROR  ANALYSIS  (GENERALIZED  MODEL) 


A.  ERROR  BOUND 


The  methodology  used  In  this  error  analysis  parallels  that  used 
In  the  development  of  equation  (III-37)  in  Section  III.  The  local 
error  bound  associated  with  the  k^h  elemental  cell  is  generated  and 
then  applied  across  the  pulse  to  develop  a  maximum  error  bound  on  a 
global  basis. 

For  a  pulse  with  arbitrary  amplitude  and  frequency  modulation,  the 
difference  between  the  actual  pulse  waveform,  s(t),  and  its  approximated 
form,  3(t)  is  given  by: 

(v-i)  s(t) -fit)  *  Wt)cosl>te>3  -vft)cos [#t>] 

where  *t)  and  ^(t)  represent  the  exact  and  approximated  instantaneous 
phase  functions  respectively  and  V(t)  and  v(t)  represent  the  exact  and 
approximated  pulse  envelopes  respectively.  The  pulse  envelope  expressions 
are  related  by: 

(V-2)  vft)  ?  v(4)  +Ve0fc) 


where  V££t)  is  the  error  that  results  from  linearly  approximating  the 
amplitude  modulation.  Making  use  of  this  relationship,  the  absolute 


pulse  error  In  the  time  domain  can  be  bounded  from  above  such  that: 


(v-3)  Ism  -  s(i)j  -  lvU)llcosl<fiuA  ^IVe/t)ljcos[^m] I 

In  the  first  term,  the  results  of  Section  (III-A)  can  be  used  to  bound 
the  absolute  difference  of  cosines  by  the  absolute  difference  of  their 
arguments.  Equations  (III-15)  thru  (III-18)  are  germane.  In  the 
second  term,  use  can  be  made  of  the  fact  that  jcofltijH  »  so  that: 

(v-4)  (sfi)-tM |  * 

Since  the  intent  is  to  extrapolate  a  local  error  to  a  global  error, 
consider  the  kth  elemental  cell  starting  at  and  of  duration,^.  If 
\Wt)  -Mdl*  ia  the  instantaneous  absolute  phase  error  associated  with 
the  kth  cell  and  Vek(t)  is  the  error  associated  with  the  amplitude 
modulation  in  the  kth  cell,  the  absolute  pulse  error  associated  with 
the  kth  cell,  |j(t)-3(t)l  can  be  written  using  equation  (IV-6) 

(v-5)  IsU)  -sit)  I*  -  l\+ ^(t  -T^) )  |  fat)  -  frt)!*  +1x^(0] 

where  A^  is  the  initial  pulse  envelope  amplitude  at  T0^  and  h*  is  the 
slope  of  the  pulse  envelope  across  the  cell.  In  a  manner  analogous  to 
equation  (III-21),  the  local  spectrum  error  associated  with  the  kth 
cell  can  be  bound  as  follows: 

(,_6)  ls(f)-M , 

If  it  is  assumed  that  the  amplitude  modulation  function,  v(t). 
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i'/K  *  4 


is  at  least  twice  differentiable,  then  in  the  elemental  cell,  the 

linear  approximation,  $(t),  results  in  error,  t),  given  by: 

(V-7)  \/tJt)  -  j  (t  -  rtK  - 

where  Tok4  V  Tok+fK  and  £  t  4  TQk+^k*  By  bounding  the  second 
derivative  of  the  modulation  function  over  the  kc^  elemental  cell, 

(V-8)  \v.Jt)\  4 

where  ■  max  ,  Tok±  tV-  •  Applying  this  bound  to 

equation  (V-6) ,  using  the  expression  for  the  absolute  instantaneous 
phase  error  developed  in  equation  (III-34),  and  replacing  the  variable 
t-Tok  by  w  the  spectrum  error  associated  with  the  kth  elemental  cell 
can  be  bounded  by: 

(v-9)  j SfV-SCf)]"  -  fj -*-MK|«2w,-3lJtwt||dw 

o  f 

t&J  |w(W-t*)Uv- 

This  expression  can  be  integrated  to  yield  the  local  bound  on  the  spectrum 
error  of  the  kc^  cell. 


(V-10) 


I 


Summing  Che  local  error  bounds  over  all  cells  provides  Che  global  specCrum 
error  bound: 


(V-ll) 


\sif)-$(f)\  *  |1  ( 

.» fr  Mcflt 


This  is  a  very  general  bound  on  Che  speccrum  error,  applicable  Co  ele- 
mental  cells  of  uniform  and  nonuniform  lengCh,  and,  as  such,  is  unwieldy. 

To  idenCify  Che  characCerisCics  of  Che  convergence  process,  some  simpll- 
ficacions  can  be  made  resulcing  in  increased  clarlcy  aC  Che  cose  of 
broadening  Che  bound.  The  expression  (MWa))  is  Che  pulse 

amplicude  aC  Che  midpoinC  of  Che  cell.  The  expression  (A^+O^h  k) 
is  Che  pulse  amplicude  aC  a  poinC  seven  eenChs  Che  widch  of  Che  kth  cell. 
BoCh  of  chese  expressions  can  be  upper  bounded  by  Che  maximum  pulse 
amplicude,  V.  Addicionally,  che  second  derivacives  of  Che  modulacing 
funcCions  can  be  bounded  by  Cheir  maximum  value  over  Che  enCire  pulse,  such 


Chac: 

(V-12a) 

/*)  s  mrfX  (^7 

■  >* 

(V-12b) 

A/  -  r*>**  Ity] 

As  a  final  simplif icadon,  lec  che  elemenCal  cells  be  of  equal  lengCh, 
f,  where  che  CoCal  pulse  widch  is  .  Wich  chese  conslderaCions, 
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equation  (V-ll)  reduces  to: 


(V-13) 


Realizing  that  2j£l  (k-1)  ■  n^-n,  the  spectral  error  bound  becomes: 
(V-14) 

Thus  the  generalized  elemental  cell  model  converges  inversely  as  the 
square  of  the  number  of  cells  used. 


B.  EMPIRICAL  ASSESSMENT 

In  order  to  assess  the  peri  .*nce  of  the  generalized  element 
cell  model  in  the  light  of  the  theoretical  error  bounds  developed, 
consider  a  pulse  with  sinusoidal  amplitude  and  frequency  modulation 
as  in  Section  (IV-B) .  The  exact  spectrum  of  this  pulse  is  known 
(Appendix  B) .  Using  the  techniques  of  Section  IV,  the  puJ.se  can  be 

modeled  using  n  elemental  cells  where  n  -  1,2,3 .  to  produce 

approximate  spectra.  For  each  choice  of  n,  the  approximate  and  exact 
spectrum  can  be  compared  to  determine  the  maximum  calculated  spectrum 
error  experienced  using  the  elemental  cell  model. 

Since  the  Instantaneous  frequency  function,  f^(t)  and  the  amplitude 
modulation  function,  v(t) ,  are  known  analytically  (Appendix  B) ,  equations 
(V-ll)  and  (V-14)  can  be  evaluated  to  generate  maximum  theoretical  error 


bounds  as  a  function  of  the  number  of  elemental  cells,  n,  used  in 
the  modeling  process.  Thus  a  set  of  three  curves  can  be  generated: 

(1)  The  maximum  theoretical  spectrum  error  bound  as  a 
function  of  the  number  of  elemental  cells  with  modulation 
variables  bounded  across  the  pulse  (equation  V-14) 

(2)  The  maximum  theoretical  spectrum  error  bound  as  a 
function  of  the  number  of  elemental  cells  with 
modulation  variables  bounded  across  elemental  cells 
(equation  V-ll) 

(3)  The  maximum  calculated  spectrum  error  as  a  function  of 
the  number  of  elemental  cells. 

The  three  curves  are  presented  in  Figure  20.  The  solid  curve 
(equation  V-14)  and  the  dashed  curve  (equation  V-ll)  display  the 
quantitative  tightening  of  the  error  bound  achieve  through  the 
increased  complexity  of  equation  (V-ll) .  The  dotted  curve  represents 
the  calculated  maximum  error  and  is  well  within  the  established 
theoretical  bounds.  In  general,  all  curves  follow  the  characteristic 
one  over  n-squared  convergence. 
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Figure  20  -  Theoretical  Error  Bound  and  Maximum  Calculated 
Errors  as  a  Function  of  the  Number  of  Elemental 
Cells  Used  to  Model  a  Pulse  with  Sinusoidal 
AM  and  FM 


VI.  EXTRAPOLATION 


A.  DEVELOPMENT 

The  concept  of  extrapolation  to  the  limit  (Richardson/Romberg 
extrapolation)  is  attractive  from  the  computational  point  of  view 
since  a  relatively  simplistic  calculation  uses  two  inexact  solutions 
to  generate  a  more  exact  solution.  The  technique  involves  some  know¬ 
ledge  of  the  order  of  the  error  in  the  inexact  solutions  and  some 
associated  assumptions.  The  desire  is  to  apply  this  type  of  technique 
to  the  elemental  cell  approach. 

Consider  the  case  of  the  unit  amplitude  pulse  with  arbitrary 
frequency  modulation  and  elemental  cells  of  uniform  length.  The 
maximum  theoretical  error  bound  for  this  case  has  been  presented  as: 

(vx-D  f--/irHfiil*eg£(i'i) 

where  again  is  the  pulsewidth,  n  is  the  number  of  elemental  cells  and 
M  is  the  maximum  value  of  the  second  derivative  of  the  instantaneous 
frequency  function  over  the  pulse.  Thus  the  approximate  spectrum  approaches 
the  exact  spectrum  roughly  as  the  inverse  of  the  number  of  cells  squared. 

Let  |P^\  be  the  approximate  spectrum  produced  when  n^  cells  are 
used  to  model  a  given  pulse,  then 


(VI-2) 


(VI-3)  /  ?(.{)- P,(f)\  ~e\ 


The  inequality  can  be  removed  by  writing: 


(VI-4) 


or  by  removing  the  absolute  value  signs: 


(VI-5)  /  P(f)l  -  jp,(f  )l  --*A  0  i.'.  4  I 


where  and  are  unknown  and  are,  in  general,  functions  of  n  and 
frequency. 

Similarly  let|Pxf  be  the  approximate  spectrum  produced  when  n£  cells 
are  used ,  then 

(VI-6)  I  PC/)  I  - 1  pjf)l  s  atxe- L  i  I 

Equations  (VI-5)  and  (VI-6)  are  actually  two  equations  in  three  unknowns: 

P,  •*/  and  aTj.  .  To  attach  some  physical  significance  to  and 
they  can  be  thought  as  of  the  normalized  errors  for  the  two  cases  (n^ 
and  n£)  discussed  above  (normalized  by  the  maximum  error  bound  for  a 
specified  number  of  cells) .  If  these  normalized  errors  (or  errors 
normalized  to  the  maximum  bound)  are  assumed  to  be  equal,  equations  (VI -5) 
and  (VI-6)  form  two  equations  in  two  unknowns.  The  equating  of  these  two 
normalized  errors  in  fact  is  an  implicit  assumption  that  for  any  given 
frequency  component  for  n^  and  n^  the  convergence  of  the  approximated 
spectrum  to  the  exact  spectrum  parallels  the  convergence  profile  established 
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by  the  maximum  theoretical  error  bound.  In  consonance  with  the 
discussions  of  Section  (III-C),  it  is  expected  that  the  validity  of 
this  assumption  is  increased  as  n  becomes  larger  and  larger.  Solving 
then  the  set  of  simultaneous  equations  for  (P|  yields: 

(VI-7)  I  PI  :  ~  C,lPJ 

(VI-8)  I  pi  r  hx  •»  J  )  1  Pa. I  ~  (Ax+  ^ 

As  is  the  case  with  extrapolation  to  the  limit  techniques,  |p\  in  equation 
(VI-8)  is  not  truly  an  exact  value  since  in  reality  ol±  is  not  equal  to 

*»• 


B.  ILLUSTRATION  AND  ASSESSMENT 

In  this  section  two  detailed  illustrations  of  the  extrapolation 
technique  are  presented  along  with  a  detailed  analysis  of  the  technique 
performance.  For  the  purposes  of  demonstration,  a  unit  pulse  with 
sinusoidal  frequency  modulation  (see  Section  II-B)  is  modeled. 

Consider  the  case  in  which  four  (4)  and  then  eight  (8)  elemental 
cells  are  used  to  model  the  pulse.  The  two  approximated  spectra  that 
result  using  the  techniques  of  Section  II  are  shown  in  figures  21  and 
22.  If  the  results  of  these  two  modeling  efforts  (n^M4;n2*8)  are  applied 
to  equation  (VI-8),  the  extrapolated  spectrum  shown  in  figure  23  is 
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Figure  21  -  Amplitude  Spectrum  for  a  Unit  Pulse  with 
Sinusoidal  Frequency  Modulation  and 
Approximated  Spectrum  Using  Four(4)  Cells 
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Figure  22  -  Amplitude  Spectrum  for  a  Unit  Pulse  with 
Sinusoidal  Frequency  Modulation  and 
Approximated  Spectrum  Using  Eight (8)  Cells 
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produced.  While  a  net  improvement  in  the  extrapolated  spectrum  can  be 
seen  in  the  first  lower  sidelobe,  certain  null  values  such  as  the 
amplitude  associated  with  f*30.12  MHz  show  deterioration.  To  gain 
a  better  insight  into  the  extrapolation  technique  performance,  figure 
24  illustrates  the  amplitude  spectrum  error  associated  with  the  four 
and  eight  cell  models  and  that  associated  with  the  extrapolated  spectrum. 

In  general,  across  the  central  region  of  the  spectrum  (f=29.4  MHz  to 
f=30.6  MHz),  the  error  associated  with  the  eight  cell  model  is  a  factor 
of  four  less  than  the  corresponding  error  associated  with  the  four  cell 
model.  Additionally,  in  this  central  region,  as  the  number  of  cells 
is  increased  from  four  to  eight,  the  convergence  to  the  exact  spectrum 
is  essentially  at  a  uniform  rate  (i.e.,  the  normalized  errors  and 

are  approximately  equal) .  The  error  associated  with  the  extrapolated 
spectrum  in  this  region  lacks  the  uniformity  seen  in  the  two  modeled 
spectra,  but  in  general  provides  a  reduction  in  the  spectrum  error  by 
about  a  factor  of  four  over  the  eight  elemental  cell  model.  On  either 
side  of  this  central  region,  the  performance  of  the  extrapolation  technique 
is  relatively  poor.  In  the  range  from  29.0  MHz  to  29.4  MHz,  the  inferior 
performance  is  the  result  of  the  fact  that  the  eight  cell  and  four  cell 
models  do  not  converge  to  the  true  spectrum  at  a  uniform  rate.  In  the 
range  from  30.6  MHz  to  31.0  MHz  the  same  problem  exists;  however,  its 
cause  is  more  graphically  evident  in  the  misalignment  of  the  spectrum 
error  sidelobes  of  the  four  cell  and  eight  cell  models. 

Before  assessing  the  performance  of  the  extrapolation  technique 
further,  consider  a  second  example  in  which  eight  cells  (n^)  and 
sixteen  cells  (n2)  are  used  to  generate  an  extrapolated  spectrum.  The 
approximated  spectra  generated  by  the  eight  cell  and  the  sixteen  cell 
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Figure  24  -  Amplitude  Spectrum  Error  for  a  Unit  Pulse  with 
Sinusoidal  Frequency  Modulation  for  Elemen  il 
Cell  Modeling  with  Four  Cells  and  Eight  Cells  and 
for  the  Extrapolation  Technique 


models  aPPear  ^  f*-8ures  22  an<*  25  respectively.  The  resultant  extra¬ 
polated  spectrum  is  not  shown  in  a  presentation  comparable  to  figure 
23  since  the  difference  between  the  extrapolated  spectrum  and  that 
produced  by  the  sixteen  cell  model  is  basically  lost  in  the  resolution 
of  the  graphical  presentation.  Instead,  the  spectrum  errors  associated 
with  an  eight  cell  model,  a  sixteen  cell  model,  and  the  extrapolated 
spectrum  (n^«*8;n2=16)  are  shown  in  figure  26.  Generally,  across  the 
spectrum,  the  convergence  of  the  sixteen  cell  model  parallels  that  of 
the  eight  cell  model  and  exhibits  the  characteristic  (1/n^)  convergence. 

The  extrapolated  spectrum  parallels  the  convergence  of  the  elemental 
cell  models  and  further  reduces  the  spectrum  error  by  about  a  factor  of 
six  or  seven.  In  this  second  example,  the  use  of  an  increased  number  of 
elemental  cells  has  reduced  or  eliminated  some  of  the  problems  addressed 
in  the  previous  example  (figure  24).  In  the  frequency  range  30.6  MHz  to 
31.0  MHz,  the  alignment  of  the  spectrum  error  sidelobes  for  the  eight 
and  the  sixteen  cell  models  has  been  greatly  improved  resulting  in 
significantly  more  accurate  extrapolated  spectrum  values.  In  the 
frequency  range  below  29.4  MHz,  convergence  of  the  two  elemental  cell 
models  is  much  more  consistent  and  anomalous  results  occur  only  between 
29  MHz  and  29.2  MHz. 

In  those  cases  when  the  assumption  of  Section  VI-A  (i.e.ff,  *  ) 

holds,  the  extrapolated  spectrum  produced  by  using  equation  (VI-8)  provides 
a  substantial  increase  in  the  accuracy  of  the  approximated  spectrum.  In 
the  case  of  the  first  example  (nj-4;n2*8),  this  is  seen  in  the  factor  of 
four  reduction  in  the  spectrum  error;  in  the  second  example  (n^*8;n2*16) , 
a  factor  of  six  reduction  in  the  spectrum  error.  The  key  point  is  that 
these  significant  enhancements  in  spectrum  accuracy  are  achieved  by 
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Figure  25  -  Amplitude  Spectrum  for  a  Unit  Pulse  with 
Sinusoidal  Frequency  Modulation  and 
Approximated  Spectrum  Using  Slxteen(16)  Cells 
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virtue  of  a  trivial  calculation.  To  form  a  basis  for  comparison  and  to 
illustrate  the  power  of  these  techniques,  figure  27  presents  an  overlay 
of  the  extrapolated  spectrum  error  (n^-8;n2“16)  discussed  above  and  the 
error  associated  with  generating  the  spectrum  with  a  4K  Discrete  Fourier 
Transform  (DFT) .  While  the  nature  of  the  error  is  different,  there  are 
few  cases  in  which  the  DFT  spectrum  provides  better  accuracy. 

The  potential  power  of  the  extrapolation  technique  is  evident  from 
the  discussions  above.  The  existence  of  anomalies  addressed  in  discussing 
figures  24  and  26,  however,  needs  further  examination.  Returning  to  the 
first  example  of  this  section  (ni*4;n2*8) ,  of  the  100  frequency  components 
between  29  MHz  and  31  MHz  examined,  there  are  32  cases  in  which  the 
extrapolated  spectrum  values  are  less  accurate  than  the  corresponding 
spectrum  values  generated  using  the  eight  cell  model.  These  32  cases 
fall  into  three  distinct  categories  which  shall  be  referred  to  as  Type(l), 
Type (2) ,  and  Type(3)  errors: 

(a)  Type(l)  Error.  In  six  instances  the  spectrum  values 
generated  using  the  eight  cell  model  are  less  accurate 
than  the  corresponding  values  generated  using  the  four 
cell  model.  Obviously,  if  |  P2 1  is  less  exact  than  |PjJ  , 
the  extrapolated  value  will,  by  equation  (VI-8) ,  be  less 
exact  than  JPjl  . 

(b)  Type (2)  Error.  In  18  instances,  the  spectrum  errors 
associated  with  {P^|  and  | P2I  are  oppositely  signed.  In 
the  extrapolation  process,  this  results  in  the  extrapolated 
value  being  driven  away  from  vice  towards  the  true  spectrum 
value.  This  type  of  error  results  from  the  nature  of  the 
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Figure  27  -  Comparison  of  Amplitude  Spectrum  Error 
Associated  with  a  4K  DFT  and  the 


Extrapolated  [n^«8;n2_16]  Elemental  Cell  Model 
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spectrum  convergence  which  will  be  discussed  below. 

(c)  Type (3)  Error.  In  eight  instances,  a  poor  extrapolated 
value  results  because  of  the  difference  in  the  rate  of 
convergence  of  the  two  approximated  spectra  to  the  true 
spectrum.  Specifically,  the  value  of  tP]J  is  considerably 
different  than  |P{,  while  the  value  of  JP^  is  a  much 
better  approximation  to  |P|  .  While  IPjl  and  I  P2I  are  not 
oppositely  signed,  the  relatively  large  difference  between 
the  two  causes,  in  essence,  an  "over-extrapolation"  to  a 
less  exact  value. 

Figure  28  presents  a  graphic  illustration  of  how  the  three  error  types 
are  actually  generated. 

To  gain  an  insight  into  the  behavior  of  the  three  types  of  error 
identified,  a  large  number  of  cases  are  examined  for  various  values 
of  n^  and  n2*  As  a  general  comment,  it  can  be  stated  that  the  char¬ 
acteristics  identified  in  the  two  examples  presented  at  the  beginning  of 
this  section  hold.  When  «c,  ft  of*, ,  the  extrapolation  technique  results  in 
a  significant  Improvement  in  the  accuracy  of  the  approximated  spectrum 
and  when  problems  occur  or  anomalies  exist,  they  are  generally  limited 
to  the  higher  ordered  sidelobes  or  local  spectrum  minima. 

Table  1  is  a  matrix  of  the  error  types  encountered  for  various 
values  of  n^  and  n£  in  exercising  the  extrapolation  technique  against 
the  unit  pulse  with  sinusoidal  frequency  modulation.  Each  block  of  the 
matrix  is  associated  with  a  specific  value  of  n^  and  n£  and  is  divided 
into  quadrants  which  represent  the  number  of  Type(l)  errors,  Type (2) 
errors  and  Type(3)  errors  as  well  as  the  total  number  of  errors  encountered 
The  legend  in  the  figure  provides  a  key  for  identifying  the  various  error 
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types.  For  example,  in  the  case  discussed  above  (n^=4;n2=8),  the  appro¬ 
priate  block  of  the  matrix  identifies  six  Type(l)  errors,  18  Type (2) 
errors,  8  Type (3)  errors,  and  3?  -otal  errors.  Examination  of  the  matrix 
yields  some  insights  into  the  different  types  of  errors  and  their  behav¬ 
ioral  characteristics.  By  examining  the  Type(l)  errors,  it  can  be  seen 
that  for  a  given  n^,  as  n2  increases,  the  number  of  Type(l)  errors  goes 
to  zero.  This  occurs  since  as  n2  increases  all  values  of  P2(f)  will 
ultimately  be  more  accurate  than  the  corresponding  value  of  P^(f)  ,  thus 
eliminating  the  source  of  the  Type(l)  error.  It  should  be  noted,  however, 
that  as  n^  increases  and  the  Type(l)  errors  diminish  that  they  may  reappear 
as  Type (2)  or  Type (3)  errors.  It  is  evident  from  the  matrix  that  the  number 
of  Type (2)  errors  is  determined  by  the  selection  of  n^.  To  some  extent 
this  is  also  true  for  Type (3)  errors,  but  there  also  tends  to  be  a 
slight  decrease  in  the  number  of  Type(3)  errors  as  n2  increases  for  a 
fixed  n^.  This  effect  results  from  the  fact  that  the  extrapolated  spectrum 
is  the  weighted  sum  of  the  two  approximate  spectra,  |P^(f)|  and  j P2 ( f ) j  > 
with  the  weights  being  determined  by  the  selection  of  nj  and  n2.  As  n2 
increases,  the  weight  associated  with  | P2 Cf ) 1  [the  more  exact  of  the  two 
approximate  spectra]  increases  and  tends  to  dominate  the  weighted  sum, 
diminishing  the  effect  of  |Pi(f)J  on  the  extrapolation  process.  In 
actual  applications,  the  relationship  between  nj  and  n2  is  critically 
important.  If  these  two  values  are  too  close  together,  the  weighted  sum 
is  biased  towards  |P^(f)|,  typically  the  less  accurate  of  the  two  appro¬ 
ximate  spectra.  On  the  other  hand,  if  n2  is  very  much  larger  than  nj_, 

|P2(f)|  may  be  so  heavily  weighted  that  the  feasibility  of  using  an 
extrapolation  approach  becomes  questionable.  Figure  29  provides  a  quan¬ 
titative  insight  into  this  problem  by  presenting  the  relative  weighting 
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Figure  29  -  Illustration  of  the  Impact  that  the  Selection  of 
n^  and  n^  Has  on  the  Weighting  of  the  Two 
Approximated  Spectrum  Values  in  the  Extrapolation 
Technique 
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of  the  spectra  |P^(f)|and  J P2 (f as  a  function  of  the  ratio  of  n^  to 
nl* 

The  most  Important  outcome  of  examining  the  matrix  In  Table  1  is 
the  realization  that  to  a  large  extent  the  successful  application  of 
the  extrapolation  technique  rests  very  heavily  on  the  proper  selection 
of  a  value  for  n^.  The  selection  of  this  value  and  the  ultimate  effect 
that  it  has  on  the  extrapolation  process  is  driven  by  the  manner  in 
which  the  approximate  spectra  converge  to  the  true  spectrum  as  the 
number  of  elemental  cells  is  varied.  To  gain  an  insight  into  the  nature 
of  this  convergence  problem,  figures  30  thru  33  present  the  convergence 
profiles  of  various  frequency  components.  The  relative  spectrum  amplitude 
[ |Pn(f)|  / (P(f)|  ]  is  plotted  as  a  function  of  n,  the  number  of  elemental 
cells  used  to  generate  the  approximate  spectrum,  )Pn(f>)  ,  for  specific 
frequency  components.  Figure  30  represents  a  well-behaved  frequency 
component  that  converges  nicely  to  the  exact  spectrum  value.  Figures 
31  thru  33,  however,  display  convergence  profiles  that  may  cause  problems 
in  applying  the  extrapolation  technique.  Figure  31  illustrates  how  a 
Type(l)  error  could  occur  with  n^»4  and  n2*8.  In  figure  32,  it  is  easy 
to  see  how  a  Type (2)  error  could  be  generated  if  n^  were  less  than  eight. 
Figure  33  presents  a  case  in  which  a  Type(3)  error  could  result  if  n^ 
were  four  and  ^  were  eight.  In  general  all  frequency  components  converge 
nicely  to  the  exact  spectrum  value  if  the  number  of  elemental  cells  is 
large  enough.  The  problem  in  dealing  with  the  extrapolation  technique 
is  to  select  n^  and  02  such  that  the  initial  transients  in  the  convergence 
profiles  do  not  impact  on  the  extrapolation  technique.  Figures  30  thru 
33  represent  three  different  types  of  convergence  prof lies  encountered  in 
examining  the  unit  pulse  with  sinusoidal  frequency  modulation.  In 


Figure  32  -  Convergence  Profile  for  f«29.9MHz 
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consonance  with  the  discussions  of  the  two  examples  presented  at  the 
beginning  of  this  section,  the  convergence  profile  displayed  in  Figure 
30  is  representative  of  those  frequency  components  within  the  malnlobe 
and  major  sidelobes.  As  frequency  components  farther  above  or  below 
the  mainlobe  are  examined,  the  convergence  profile  becomes  more  complex. 

The  extrapolation  technique  has  been  observed  to  provide  significant 
increases  in  the  accuracy  of  spectrum  estimates  at  a  m-fn-taml  cost  in 
computation.  The  use  of  this  technique,  however,  requires  a  judicious 


selection  of  n-i 


VII.  APPLICATION 


The  model  developed  in  the  previous  sections  can  be  a  useful  tool  in 
the  development  and  implementation  of  signal  classification  algorithms 
based  on  frequency  domain  analysis.  The  modulation  that  is  present  on 
pulsed  signals  may  well  be  a  useful  discriminant  in  the  classification 
process.  This  being  the  case  the  treatment  of  this  discriminant  in  the 
frequency  domain  has  the  potential  of  being  less  susceptible  to  noise 
since  the  transformation  is  a  smoothing  process. 

To  address  the  discriminant  selection  in  the  frequency  domain,  con¬ 
sider  the  chirp  radar  pulse.  As  the  chirp  rate  is  increased  or  decreased 
for  a  given  pulse  duration,  certain  characteristic  effects  occur:  (I) 
the  spectrum  is  shifted  in  proportion  to  the  chirp  rate  and  (2)  the  side- 
lobe  amplitudes  change  relative  to  the  mainlobe  as  a  function  of  chirp 
rate.  Knowledge  of  the  form  of  the  spectrum  provides  an  Insight  into  the 
nature  of  the  chirp  pulse.  In  more  general  cases  in  which  the  frequency 
modulation  may  be  more  complex,  does  knowledge  of  the  spectrum  form  pro¬ 
vide  useful  information  to  the  classification  process? 

To  examine  this  question,  let  the  discriminants  that  define  the 
spectrum  form  do  so  in  terms  of  the  sidelobe  structure  relative  to  the 
mainlobe.  This  can  be  easily  accomplished  by  choosing  two  types  of 
discriminants:  (1)  the  relative  amplitudes  of  the  sidelobe  local  maxima 
(relative  to  the  maximum  spectrum  amplitude)  and  (2)  the  frequency  dif¬ 
ferences  between  the  sidelobe  maxima  and  the  center  frequency  of  the 
mainlobe.  The  selection  of  these  parameters  establish  the  form  of  the 
spectrum. 

To  actually  assess  the  feasibility  of  employing  sidelobe  structure 
in  the  classification  process  requires  that  the  Impact  of  the  modulation 
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on  Che  spectrum  be  understood.  To  gain  an  insight  into  this  effect,  con¬ 
sider  a  unit  pulse  with  an  instantaneous  frequency  function  as  shown  in 
figure  34.  Note  that  the  initial  portion  of  the  pulse  (0.2/cs)  is  subject 
to  a  symmetrical  triangular  modulation  of  magnitude  4FM.  To  investigate 
the  impact  of  the  modulation  on  the  frequency  domain,  the  values  of  the 
two  discriminants  (relative  amplitude  and  relative  frequency)  are  deter¬ 
mined  for  the  first  three  upper  and  lower  sidelobes  for  a  variety  of 
pulses  as  AFM  varies  from  zero  to  5  MHz.  Figure  35  displays  the  relative 
amplitudes  of  the  upper  sidelobe  maxima  as  a  function  of  AFM.  For  exam¬ 
ple,  the  first  sidelobe  has  a  relative  maxima  of  0.3  when  the  magnitude 
of  the  triangular  modulation  is  aFM  ■  2.0  MHz.  All  three  sidelobe  maxima 
have  the  potential  of  functioning  as  discriminants  since  they  are  single 
valued  and  monotonically  Increasing.  For  example,  the  relative  amplitude 
of  the  sidelobe  for  AFM  -  2.0  MHz  is  easily  distinguished  from  that  for 
AFM  *2.5  MHz  for  all  three  sidelobes.  Figure  36  is  a  similar  presenta¬ 
tion  for  the  lower  sidelobe  maxima.  For  a  pulse  with  initial  triangular 
modulation,  the  lower  sidelobes  are  inferior  to  the  upper  sidelobes  as 
discriminants.  This  is  due  to  the  fact  that  the  relative  amplitudes 
extend  over  a  much  more  confined  range  of  values  and  that,  in  the  case  of 
the  second  sidelobe,  points  exist  on  the  curve  that  are  multi-valued. 

From  figure  36  it  appears  that  the  first  sidelobe  disappears  when  AFM  - 
2.5  MHz.  In  fact,  as  the  modulation  magnitude  increases  past  2.0  MHz, 
the  first  and  second  sidelobes  merge  to  form  one  and  are  arbitrarily 
associated  with  the  second  sidelobe.  The  inferior  performance  of  the 
lower  sidelobes  in  this  example  should  not  be  cause  to  disregard  them  as 
discriminants,  since  by  changing  the  polarity  and  location  of  the  trian¬ 
gular  modulation  the  relative  merits  of  the  upper  and  lower  sidelobes 
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Figure  34  -  Instantaneous  Frequency  Function  with 
Triangular  Modulation 
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Figure  35  -  Upper  Sidelobe  Relative  Amplitude  vs. 
Frequency  Modulation 
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Relative  Amplitude  of 
0.5  f  Lower  Sidelobe  Maxima 


(a)  First  Sidelobe 

(b)  Second  Sidelobe 

(c)  Third  Sidelobe 
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Figure  36  -  Lower  Sidelobe  Relative  Amplitude  vs. 
Frequency  Modulation 


would  be  reversed.  The  relative  amplitude  of  the  sidelobe  maxima  display 
the  characteristics  desired  of  a  discriminant  for  classification  pur¬ 
poses. 

Figure  37  displays  the  relative  frequency  of  the  upper  sidelobe 
maxima  as  a  function  of  4FM.  In  order  to  accomodate  all  three  sidelobe 
profiles  on  the  same  scale,  fQ  is  selected  as  the  normalization  factor. 

It  represents  the  nominal  frequency  at  which  the  sidelobe  maxima  would 
occur  for  the  case  of  a  constant  frequency  pulse.  If  f  is  the  pulse 
duration,  the  sin(x)/x  spectrum  will  have  local  maxima  at  the  following 
approximate  spacings  from  the  center  of  the  mainlobe:  +3/2 f,  +5/2f,  +7/2^*, 
...  so  that  f0,  the  location  of  the  nCh  sidelobe  peak  relative  to 
the  center  frequency,  is  +(n  +  l/2)/f.  Thus  as  each  of  the  three  side- 
lobes  is  considered,  the  value  of  fQ  changes.  The  range  of  values  over 
which  the  relative  frequencies  extend  is  fairly  small  especially  for  4FM 
less  than  2.0  MHz.  Figure  38  presents  the  relative  frequency  of  the 
lower  sidelobe  maxima.  The  first  sidelobe  provides  no  useful  informa¬ 
tion.  The  second  and  third  sidelobes  however  are  quite  robust  due  in 
large  part  to  the  fact  that  the  first  and  second  sidelobes  merge  at  about 
aFM  ■  2.0  MHz.  There  appears  some  potential  for  relative  frequency  as  a 
discriminant  but  its  power  does  not  appear  comparable  to  that  of  relative 
amplitude. 

To  gain  another  feel  for  discriminant  power  consider  a  pulse  as  in 
figure  3A  but  with  a  triangular  modulation  that  is  not  symmetrical.  Such 
an  instantaneous  frequency  profile  is  shown  in  figure  39.  Again,  the 
relative  amplitude  and  relative  frequency  of  the  local  sidelobe  maxima 
can  be  examined.  Figure  AO  displays  the  relative  amplitude  for  the  upper 
sidelobe  maxima.  The  results  from  figure  35  (symmetrical  modulation)  are 
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Figure  37  -  Upper  Sidelobe  Relative  Frequency  vs. 
Frequency  Modulation 
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Figure  38  -  Lower  Sidelobe  Relative  Frequency  vs.  Frequency 
Modulation 
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plotted  as  solid  lines  for  comparison.  These  discriminants  appear  to 
work  well  in  distinguishing  between  these  similar  pulses  especially  for 
AFM  greater  than  2.5  MHz.  A  similar  presentation  for  the  lower  sidelobe 
relative  amplitudes  appears  in  figure  41.  The  performance  of  the  lower 
sidelobe  discriminants  parallels  that  of  the  upper  sidelobes  for  these 
two  similar  pulses.  Figures  42  and  43  present  the  relative  frequency  of 
the  upper  and  lower  sidelobes,  respectively,  as  a  function  of  4FM. 

Except  for  large  values  of  dFM  in  the  case  of  the  second  and  third  lower 
sidelobes,  the  overall  performance  of  relative  frequency  in  distin¬ 
guishing  between  these  two  signals  is  marginal. 

In  summary,  frequency  domain  techniques  are  reasonable  for  pulse 
classification  purposes.  The  use  of  sidelobe  structure  in  the  form  of 
relative  amplitude  and  relative  frequency  as  a  classification  discrimi¬ 
nant  is  a  valid  approach  to  the  problem  with  the  former  possessing  the 
greater  potential  for  success. 

This  type  of  analysis  lends  itself  very  well  to  the  elemental  cell 
model.  Waveforms  are  easily  constructed  and  implemented  and  local  maxima 
can  be  readily  identified  to  arbitrary  precision. 


112 


Relative  Amplitude  of 
Lower  Sidelobe  Maxima 


Symmetrical  Modulation 
Asymmetrical  Modulation 


4L 


(a)  First  Sidelobe 
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Figure  42  -  Upper  Sldelobe  Relative  Frequency  vs. 
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Figure  43  -  Lower  Sidelobe  Relative  Frequency  vs. 


Frequency  Modulation 


VIII.  CONCLUSIONS  AND  RECOMMENDATIONS 


A  flexible,  easy  Co  use  tool  for  modeling  the  spectra  of  pulse 
signals  with  arbitrary  amplitude  and  frequency  modulation  has  been 
developed  along  with  the  appropriate  error  bounds.  When  properly 
^implemented,  accuracy  comparable  to  that  of  an  FFT  is  achievable. 

The  technique  developed  is  not  universally  superior  to  the  FFT; 

efficiency  of  one  over  the  other  is  a  function  of  specific  applications. 

Computational  spded  can  be  an  advantage  in  the  elemental  cell  method 
in  those  cases  where  large  FFT's  are  required  relative  to  the  product 

of  the  number  of  elemental  cells  (n)  and  the  number  of  frequency 

components  required  (m) . 

The  feasibility  of  using  frequency  domain  discriminants  in  the 
signal  classification  process  has  been  demonstrated. 

There  are  a  number  of  areas  for  further  investigation  and  refinement. 
Appendix  D  addresses  the  comparison  of  the  FFT  and  the  elemental  cell 
technique  in  terms  of  computational  speed.  Since  most  of  the  studies 
•  during  this  effort  were  focused  on  precision  vice  speed  and  since  the 
execution  time  studies  were  conducted  without  modification  to  the  basic 
subroutines  (Appendix  A),  significant  Improvements  can  be  made  in  the 
efficiency  of  the  elemental  cell  subroutines.  Several  specific  areas 
for  improvement  are: 

(1)  the  conversion  of  double  precision  based  routines  to 
single  precision.  Preliminary  investigation  in  this  area  has  indicated 
that  this  can  be  accomplished  with  little  degradation  in  model 
performance. 
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(2)  the  elimination  of  “left-hand"  spectrum  computations.  The 
existing  subroutines  calculate  both  the  “left-hand"  and  “right-hand" 
spectrum  contributions  to  each  frequency  component.  Under  the  proper 
conditions,  the  former  can  be  eliminated  without  performance  degradation. 

(3)  the  streamlining  of  certain  aspects  of  the  overall  computation. 
For  example,  the  technique  of  James  (5),  discussed  in  Appendix  C,  for 
Fresnel  integral  evaluation  should  be  assessed  for  use  in  place  of  the 
existing  technique. 

The  time  required  to  generate  a  spectrum  using  the  elemental  cell 
technique  is  a  function  of  the  number  of  cells  used  (n)  and  the  number  of 
frequency  components  required  (m).  The  time  required  increases  as  the 
product  of  m  and  n.  Since  the  calculation  of  each  elemental  cell  spectra 
is  Independent  of  all  other  elemental  cell  spectra,  individual  elemental 
cell  spectra  can  be  calculated  in  parallel  and  then  coherently  summed. 
Investigation  of  the  implementation  of  parallel  processing  arrays  should 
be  initiated.  Under  such  an  implementation,  the  computation  time 
required  for  spectrum  generation  would  be  proportional  to  m  vice  m  x  n. 

While  the  feasibility  of  using  spectrum  discriminants  in  signal 
classification  has  been  established,  further  examination  is  required: 

(1)  the  approach  should  be  assessed  against  a  noisy  environment; 

(2)  the  approach  should  be  compared  to  other  techniques  on  simulated 

and  real  data. 

Further  work  might  well  be  done  in  applying  the  elemental  cell 

concept  with  different  basis  functions  that  lend  themselves  to  certain 

classes  of  signals. 
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APPENDIX  B 


SPECTRA  OF  SINUSOIDALLY  MODULATED  PULSES 


Consider  a  pulse,  s(t),  which  may  be  both  amplitude  and  frequency 
modulated  so  that  over  the  duration  of  the  pulse: 

(B-l)  sfc)  =  v/(t)cos 

where  v(t)  represents  the  amplitude  modulation  and  *^(t) ,  the  instanta¬ 
neous  phase,  accounts  for  the  frequency  modulation.  An  expression  for 
the  frequency  spectra  of  two  specific  cases  of  equation  (B-l)  is 
developed.  In  CASE  A,  the  pulse  envelope  is  simply  represented  by  a 
unit  pulse  (i.e.,  v(t)-l  over  the  pulse  width).  In  CASE  B,  the  pulse 
envelope  is  allowed  to  vary  sinusoidally.  In  both  CASES  A  and  B,  however 
the  carrier  frequency  experiences  a  sinusoidal  frequency  modulation. 

Consider  a  sinusoidal  instantaneous  frequency  function  as  shown  in 
figure  (B-l)  with  fc  being  the  carrier  frequency;  ^m*  the  modulating 
frequency  and  kf ,  the  modulation  index.  The  instantaneous  frequency 
function,  f^(t),  is  represented  by: 

(B-2)  f;(t)  --fc  t 

The  instantaneous  phr.se  can  then  be  written  as: 


(B-3) 


^Kt)  =  (JJ9t  t 


where  ^(o)  is  the  initial  carrier  phase  at  time  zero.  For  the  purposes 
of  this  analysis,  the  interest  is  in  the  cosine  of  the  phase  function 

ft  t) 

(B-4)  Cosl^til  -  COs[wet  *  &fS/AJ  (oJfint)  f  ^Coij 


The  right  hand  side  of  equation  (B-4)  can  be  expanded  as  an  infinite 
series  of  sinusoidal  functions  weighted  by  Bessel  functions  (6) . 


(B-5) 


cos 

+(-l)1  COS 


This  expression  represents  the  sinusoidal  FM  associated  with  the  pulse 
for  which  spectra  are  to  be  generated  for  CASES  A  and  B  below. 


CASE  A:  SINUSOIDAL  FM 


The  frequency  spectrum  for  the  unit  pulse  of  duration^ith  sinusoidal 
FM  can  be  expressed  as: 

ZCf)  *  f  cos 


(B-6) 


Using  the  results  of  equation  (B-S)  and  integrating,  produces: 


(B-7) 


Z(f)  -  /]- 


-J&o) 


f gjCiMc-iui/m*'*)? 


*(z£k^l  ] 

uju  r  iu>/tm  *  60  I?  J 


i-iUJ/vm.  i 


CUe/)tU/»M 

This  frequency  representation  is  essentially  the  convolution  of  the 
spectrum  of  a  sinusoidally  modulated  carrier  with  the  spectrum  of  a 
unit  pulse.  That  is  a  sine(x)/x  spectrum  convolved  with  an  infinite 
sequence  of  delta  functions  weighted  by  the  Bessel  functions,  J^kf). 

Thus  the  spectrum  is  a  function  of  the  modulation  index  k^.  If  k^  is 

small,  the  first  few  delta  functions  account  for  most  of  the  power  in 

\ 

the  carrier(6).  For  example,  if  kj»l,  then  Jq,  J^,  and  J2  account  for 
98%  of  the  signal  power.  In  using  equation  (B-7)  to  generate  a  spectrum 
for  CASE  A,  the  first  twenty  terms  of  the  series  are  evaluated.  A  program 
listing  which  evaluates  equation  (B-7)  appears  at  the  rear  of  the  Appendix 
and  is  titled  "SUBROUTINE  SINE". 


B.  CASE  B:  SINUSOIDAL  AM  AND  FM 


Consider  now  CASE  B  in  which  Che  amplitude  modulation  is  allowed  to 
assume  a  sinusoidal  form  as  shown  in  figure  B-2.  The  frequency 
modulation  is  unchanged  and  generates  the  cos  0(t)  term  listed  in 
equation  (B-5).  Allow  the  amplitude  modulation  to  be  expressed  as: 


(B-8) 


V/(t)  =  S/Jd 


where  f  *  l/l'f  .  Thus  the  pulse  waveform  equation  can  now  be  written: 

cl 


(B-9) 


s/i)-  S/A/ 


J  Coi  [jpte) 3 


With  the  aid  of  equation  i.B-4)  and  trigonometric  identities,  the  pulse 
equation  can  oe  expressed  as: 


(B-10) 


s(t)  -  ^  ^ ^3  *  S,AJ 

■f  21  j  S //0[(tue  -  i  oj/m  * 

-Si a/ [fouc ~ luJ/m 

+ (- f)ls/io  fo*  +  i  )d 

+  S/A/[(uic+i  «/** 

By  taking  the  Fourier  transform  of  the  pulse,  the  frequency  spectrum 
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for  CASE  B  can  be  expressed  as: 


f  c-jW# {  f  [, 

\.  "«- *«'/•«  ♦i'O'Ui.***'  J 

*(  if*'*1  4’('0<UJ^u^t_  |\ 

\  Wfc *tu»/>n4  C-Oco^-U*  J 


K+i  t  /  /  '<»»*  *  C’<)W  +ujf  j 

\  UJe*l  CO/yt,  +(-!./  U>«, 4-  U>  „ 


Again  as  with  CASE  A,  the  first  twenty  terms  of  the  series  are  evaluated 
In  establishing  a  reference.  A  program  listing  for  CASE  B  appears  at 
the  rear  of  the  appendix  and  is  titled  "SUBROUTINE  ASINE". 
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CALCULATE  THE  FIRST  TWENTY  TERMS  OF  THE  BESSEL  SERIES  EXPANSION 
AND  FORM  THE  CUMULATIVE  SUM 
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APPENDIX  C 


FRESNEL  INTEGRALS 

The  Fresnel  integral,  important  in  the  theory  of  diffraction,  can 
be  written  as: 

(c-i)  C(*)-jScx)  -  J* 

Plotting  this  integral  in  the  complex  plane  yields  a  curve  known  as 
Cornu's  spiral  (figure  C-l) .  The  positive  values  of  x  appear  in  the 
first  quadrant  and  the  negative  values  in  the  third.  A  vector  drawn  from 
the  origin  to  any  point  on  the  curve  represents  in  both  phase  and 
magnitude  the  value  of  the  integral  in  equation  (C-l) .  (The  phase  of  the 
vector  is  the  negative  of  the  phase  of  the  integral.)  Note  from  the 
Cornu  spiral  that  as  x  goes  to  infinity,  C(x)  and  S(x)  both  converge 
to  one  half. (7) 

The  integral  in  equation  (C-l)  can  be  expressed  in  terms  of  a  series 
expansion  as: 

(C-2a)  CCt)  =  2-  (lH)!(Wn+0 

(C-2b)  S(t)-  2-  QmQHto+J) 

The  fact  that  the  Fresnel  integral  is  not  expressible  in  closdd  form 
means  that  in  application  approximations  to  equation  (C-l)  must  be 
employed.  One  such  approach  has  been  developed  by  Boersma(2) 
in  which  equation  (C-l)  is  written  with  a  change  of  variable  (w-^v^) 


(C-3) 


Cu)-yJa)  *  [* 

1  Jo 


Boersma's  technique  makes  use  of  the  following  approximations: 

(c-4a>  Ca)-jS(i')  -  e**  £/*vW@ 

_ _  n  w 

(c-4b)  C(i)-jSa)--  +  Z.(c„>jd*)ii)  2* 

'  tuO  ° 

where  the  coefficients  afl,  bQ,  cQ  and  dn  are  derived  by  the  f-method 
of  Lanczos  and  are  provided  in  (2).  The  maximum  error  associated  with 
this  approach  is  1.6x10“^.  This  technique  of  evaluating  the  Fresnel 
Integrals  has  been  employed  in  the  subroutine  CS  listed  in  Appendix  A. 

Other  techniques  and  approximations  are  also  available  to  evaluate 
the  Fresnel  integrals;  one  of  the  more  recent  being  presented  by 
James (5).  James'  approach  is  unique  in  that  it  employs  one  concise 
approximate  expression  for  all  real  values  of  the  argument.  The 
technique  is  implementable  on  a  pocket  calculator  and  provides  a 
maximum  amplitude  error  of  eight  per  cent. 


APPENDIX  D 


COMPARISON  OF  FFT  AND  ELEMENTAL  CELL  EFFICIENCY 

Given  a  sequence  x(n)  where  n  ■  1,2 . (N-l),  the  Fourier 

coefficients  associated  with  the  sequence  can  be  represented  by  the 
discrete  Fourier  transform  (DFT): 

(D-i)  Xfi)  =  X(n)  VJfii  o±A±N-I 

nto 

.;«£) 

where  Wn  .  A  straightforward  calculation  using  equation  (D-l) 

would  require  N2  operations  where  the  term  "operation"  is  taken  to 
mean  a  complex  multiplication  followed  by  a  complex  addition.  More 
efficient  methods  are  available  for  evaluating  equation  (D-l)  which 
require  considerably  fewer  than  N2  operations.  Cooley  and  Tukey(8) 
developed  a  technique  which  has  lead  to  a  class  of  algorithms  known 
as  Fast  Fourier  Transforms  (FFT's)  which  provided  for  the  evaluation 
of  equation  (D-l)  in  less  than  2Nlog2N  operations. 

While  the  number  of  operations  required  is  frequently  used  in  DFT 
and  FFT  analysis  as  a  measure  of  computational  complexity  or  efficiency 
it  cannot  easily  be  applied  to  the  elemental  cell  model  to  establish 
a  common  base  for  comparison.  The  actual  time  required  for  an  FFT 
or  elemental  cell  calculation  can  be  compared. 

Using  the  subroutine  HARM  (IBM  Scientific  Subroutine  Library) , 
which  is  based  on  the  algorithm  of  Cooley  and  Tukey(8),  the  Fourier 
transforms  of  50  signals  were  generated  to  determine  the  average 
execution  time  of  an  N-polnt  transform  where  N  -  512,  1024,  etc. 
for  each  set  of  50  signals. 


In  similar  fashion,  the  elemental  cell  spectra  were  generated  for 
each  of  the  four  classes  of  elemental  cells  using  one  hundred  frequency 
components.  Table  (0-1)  lists  the  results  of  this  analysis.  In  the 
case  of  the  FFT,  as  the  size  of  the  transform  is  doubled,  the 
time  required  for  computation  increases  by  approximately  a  factor  of 
two.  This  is  in  keeping  with  what  would  be  predicted  by  the  measure 
of  complexity  discussed  above.  In  the  case  of  the  elemental  cell  model 
the  time  required  for  computation  Increases  with  the  sophistication 
of  the  class  of  elemental  cell. 

In  application,  the  time  required  to  implement  an  FFT  is  a  function 
of  the  transform  size(N),  which  is  in  turn  a  function  of  the  problem 
being  addressed.  The  time  required  to  implement  the  elemental  cell- 
model,  on  the  other  hand,  is  not  only  a  function  of  the  problem 
being  addressed  but  also  the  nature  of  the  elemental  cells  employed, 
the  number  of  cells  used  and  the  number  of  frequency  components  to 
be  generated.  For  example,  the  average  execution  time  for  an 
elemental  cell  model  composed  of  ten  constant  amplitude  cells  with 
linear  FM  for  which  twenty  spectral  components  are  desired  is 
(10  x  20  x  6.10msec  «)  1.22  seconds. 

A  graphical  comparison  of  the  execution  times  of  the  two 
techniques  appears  in  figure  (D-l) ,  in  which  execution  time  is 
plotted  as  a  function  of  the  product  of  the  number  of  elemental 
cells  (n)  and  the  number  of  frequency  components  required  (m) . 

Nothing  can  be  said  unequivocally  about  the  computational  efficiency 
of  the  elemental  cell  model  vis-a-vis  the  FFT  apptoach.  Under 
certain  circumstances,  when  the  transform  size  is  large  (e.g.  greater 
than  2K)  and  the  product  of  m  and  n  is  small  (e.g.  less  than  200) 


AVERAGE  TIME  TO  EVALUATE  AN  N -POINT  FFT 


_N _ 

TIME 

512 

0.22 

sec 

1024 

0.47 

sec 

2048 

0.88 

sec 

4096 

1.81 

sec 

8192 

4.01 

sec 

AVERAGE  TIME  TO  EVALUATE  AN  ELEMENTAL  CELL  FOR  ONE  FREQUENCY 


CLASS  OF  CELL 


CONSTANT  AM 
LINEAR  AM 
CONSTANT  AM 
LINEAR  AM 


CONSTANT  IM 
CONSTANT  FM 
LINEAR  FM 
LINEAR  FM 


TIME 

3.72  msec 
5.55  msec 
6.10  msec 
8.85  msec 


Table  D-l  -  Average  Times  for  FFT  and  Elemental  Cell 


Calculations 


n  ■  number  of  cells 

m  ■  number  of  frequency  components  desired 
Figure  D-l  -  FFT  and  Elemental  Cell  Execution  Times 
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speed  can  be  obtained  using  the  cell  model  at  the  cost  of  discriminating 
the  signal  of  interest. 

It  should  be  noted  that  in  this  analysis  the  subroutines  used  to 
generate  the  elemental  cell  spectra  are  those  listed  in  Appendix  A 
and  that  no  attempt  at  time  optimization  has  been  made.  In  fact, 
at  the  cost  of  efficiency,  precision  has  been  stressed.  It  is  expected 
that  by  modifying  the  existing  subroutines  to  handle  only  the  "right 
hand"  spectrum  and  to  function  in  single  vice  double  precision  that  the 
execution  time  could  be  reduced  by  nearly  a  factor  of  four. 
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